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NOTE. 



The following pages are printed from the stereotype plates 
of the Appendix to the author's Manvxil of Spherical and 
Praeiical Astronomy^ without any change in the numbering 
of the pages or of the tables. The text, therefore, com- 
mences with p. 469 and ends with p. 566 ; and the tables 
are numbered ix., ix. a., x. and x. A., as they stand in the 
original work. 

As the Method of Least Squares is applicable in almost all 
the physical sciences where numerical results are to be de- 
daced Lm observations, and i, here t«.ted ftom fund- 
mental and general principles, it is thought that this separate 
publication of the Appendix will supply the want that has 
for a long time been felt of a suitable text-book on this sub- 
ject for the use of students of practical science generally, and 
more especially of classes in our scientific schools. 

St. Louis, January 1, 1868. 




METHOD OF LEAST SQFARE3.* 



1. A NUMBER of obflerratioiia being taken for the purpose of 
determining one or more unknown quantities, and these obser- 
vations giving discordant resulta, it is an important problem to 
determine the most probable values of the unknown quantities. 
The method of least squares may be defined to be that method 
of treating this general problem which takes as its fundamental 
principle, that the most probable values are those wkkk make the sum 
of the squares of the residual errors a minimum. But, to understand 
this definition, some degree of acquaintance with the method 
itaelf ia necessary. 

* The tret publiehed application of the method is to be found in Leqendbe, Nouvrlla 
mfthadei pour la Mtermmation da orbites dts comitri, Paris, 1806. The development, 
hoTever, ft'otn fuudamenta,! prinaiplea is du« to OAuaa, who declared that he bad 
used the method as eart; ub 1TB5. See his Tktaria Motut Corpomm Caliatium, 1809, 
Lib. n. Sbo. III.; Diiquailio de elemmtia ell^Udi Palladii, 1811; Beslimmmg der 
Oenanigkcil der Beohachtmigen (t. Lindenav and BoHNENBKsoBB'a ZaltchTift, 1816, 1. 
B. 16G) i Theoria comhinationit obfervalionttm erronbut mimmis abnoxiK, 1823 ; Supple- 
KenUim theeria coTnbinationU, &c,, 1826: ^U of which have been rendered qnite acceae- 
iblo through a. Frenoh tronslalion by J. Bertband, Melhodt dts moindret carriii. M(- 
fnsirei mr la combination det obiervations, par Cs. Fr. Gaubh, Paris, 1855. 

For a digest of the preceding, together with Che reanlta of the labors of Bessel 
and HakbbN, see Encke, I'eber die Melhade der Heimlen Quadrale, Berliner Aatron. 
Jahrbuch for 1834, 1835, 183S: in oonneotioD with which must be mentionpd cape- 
oiiilly the practical worlt of GBBllNa, -Die AuBgltiehungtreeknungen der praetiteken 
Qeomeirit, Hamburg, 1843, 

Bee also Iaplaob, Thiorie rmxlgtiqui dew probabmift, LiT. 11. Chap. IV. ; Poisson, 
5ur laprobabilii/des rttultali moijem da obtervaliam, in (he CannaissaDce dea Temps for 
1827 ; Encke, in the Berlin Jahrbuoh for 1853 ; Besbei., In Attron. Nach., Xos. 358, 
859. 399; Hansen, in AiCron. Nofh., Nob. 192, 202 et seq. i Pbiroe, in the Attron, 
Jo-amal (Caoihridge, Mass.), Vol. 11. No. 21 ; Xiaobe, ChIchI da probMlitii tt tMorit 
da tTTnin, Bruxellea, 1862. 
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ERRORS TO WHICH OBSBBVATIONS ARK LIABLE. 



. Kvery obBerration which is a measure, however carefully i 
may be made, is to be regarded as subject to error; for expe 
rieiice teaches that repeated measurea of the same quantity, loAefl 
the greatest precision is sought,* do not give uniformly the eami 
result. Two kinds of errors are to be diBtinguished. 

Chnatant or regular errors are those which in all measures of th< 
same quantity, made under the same circumstances, obtain th( 
same magnitude; or whose magnitude is dependent upon the 
circumstances according to any determinate law. The causes ol"* 
such errors must be the subject of careful preliminary search m 
all physical inquiries, so that their action may be altogether pre- 
vented or their effect removed by calculation. For examplfl^y 
among the constant errors may be enumerated refraction, aber- 
ration, &.C. ; the effect of the temperature of rods used in i 
Buring a base line in a survey ; the error of division of a graduated 
instrument when the same division is used in all the measurea ; 
any peculiarity of an instrument which affects a particular mea> 
eurement always by the same amount, such as inequality of tho 
pivots of a transit instrument, defective atljustment of the colli- 
mation, imperfections of lenses, defects of micrometer screws, &e. 
to which must be added constant peculiarities of the observer, 
who, for example, may always note the passage of a star over a. 
thread of a transit instrument too soon, or too late, by a eonstani 
quantity, or who, in attempting to bisect a star with a mierometM 
thread, constantly makes the upper or the lower portion th( 
greater; or who, in observing the contact of two images (in 
sextant measures, for instance), assumes for a contact a position 
in which the images are really at some constant small distaucsj 
or a position in which the images are really overlapped, &e. &o, 

Thus, we have three kinds of constant errors : 

Ist. Theorelieal, such as refraction, aberration, &c., whose effects, , 
when their causes are once thoroughly understood, may be cal- 
culated a priori, and which thenceforth cease to exist as errors. 

* The qualification, " when the greatest preoision is sought," ie important ; for if, 
e.j., we were to determine the latitude of a place liy repeated menBures of the raeri- 
diau altitude of the same fixed stnr with a eextatit diTided onl; to whole degrees, all 
our measareB might give the snme degree. The nacordunce of obaervations is, there- 
fore, not to be taken as an infallible evidence of tlieir nccuraoj. It ib especially 
when we approach t/ie limits of our measurinff paacri that we beoouie sensible o 
dixcrepuicies of observationa. 
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The detection of a constant error in a certain class of observa- 
tions very commonly leads to iuveatigatiuns by which its cause 
is revealed, and thus our physical theories are improved. 

2d. InsCntmental, which are discovered by an examination of 
our instruments, or from a diacus&ion of the observations made 
with them. These may also be removed when their causes are 
fully understood, either by a proper mode of using the instru- 
ment, or by subsequent computation. 

3d. Personal, which depend upon peculiarities of the observer, 
and in delicate inquiries become the subject of special investigar- 
tion under the name of "personal equations." 

We are to assume that, in any inquiry, all the sources of con- 
stant error have been carefully investigated, and their eftects 
eliminated as far as practicable. When this has been done, 
however, we find by experience that there still remain discrepan- 
cies, which must be referred to the next following class. 

Irreffular or accidental errors are those which have irregular 
causes, or whose effects upon individual observations are gov- 
erned by no fixed law connecting them with the circumstances 
of the observations, and, therefore, can never be subjected 
a priori to computation. Such, for example, are errors arising 
from tremors of a telescope produced by the wind ; errors in the 
refraction produced by anomalous changes of density of the 
strata of the atmospliere ; from unavoidable changes in the 
several parts of an instrument produced by anomalous variations 
of temperature, or anomalous contraction and expansion of the 
parts of an instrument even at known temperatures ; but, more 
especially, eri'ors aiiaing from the imperfection of the senses, as 
the imperfection of the eye in measuring very small spaces, of 
the ear in estimating small intervals of time, of the touch in the 
delicate handling of an instrument, &c. 

This distinction between constant and irregular errors is, 
indeed, to a certain extent, rather relative than absolute, and 
depends upon the sense, more or less restricted, in which we 
consider observations to be of the same nature or made under the 
same drcumstanees. For example, the errors of division of an 
instrument may be regarded as constant errors when the same 
division comes into all measures of the same quantity, but as 
irregular when in every measure s, different division is used, or 
when the same quantity is measured repeatedly with different 
InstrumeutB. 
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After a full investigation of the constant or regular errors, it 
is the next business of the observer to diminish as much as pos- 
sible the irregular errors by the greatest care in the observations ; 
and finally, when the observations are completed, there remains 
the important operation of combining them, so that the outstand- 
ing, unavoidable, irregular errors may have the least probable 
effect upon the results. For this combination we invoke the 
aid of the method of least squares, which may be said to have 
for its object the restriction of the effect of irregular errors within 
Ihe narrowest limits according to the theory of probabilities, and, 
at the same time, to determine from the observations themselves 
the errors to which our results are probably liable. It is proper 
to observe here, however, to guard against fallacious applications, 
that the theory of the method is grounded upon the hypothesis 
that we have taken a large number of observations, or, at least, a 
number suflBlciently large to determine the errors to which the 
observations are liable. 

CORRECTION OP THE OBSERVATIONS. 

8. When no more observations are taken than are sufficient 
to determine one value of each of the unknown quantities 
sought, we have no means of judging of the correctness of the 
results, and, in the absence of other information, are compelled 
to accept these results as true, or, at least, as the most probable. 
But when additional observations are taken, leading to different 
results, we can no longer unconditionally accept any one result 
as true, since each must be regarded as contradicting the others. 
The results cannot all be true, and are all probably, in a strict 
sense, false. The absolutely true value of the quantity sought by 
observation must, in general, be regarded as beyond our reach ; 
and instead of it we must accept a value which may or may not 
agree with any one of the observations, but which is rendered 
most probable by the existence of these observations. 

The condition under which such a probable value is to be 
determined, is that all contradiction among the observations is to be 
removed. This is a logical necessity, since we cannot accept for 
truth that which is contradictory or leads to contradictory results. 

The contradiction is obviously to be removed by applying to 
the several observations (or conceiving to be applied) probable 
eorrectionSj which shall make them agree with each other, and 
which we have reason to suppose to be equivalent in amount to 
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METHOD OF LEAST SQUARES. 4T3 

the accidental errors severally. But let ub here reniark that we 
do uot in this statement by any means imply that an observer is 
to arbitrarily assume a system of corrections which will produce 
accordance : on the contrary, the method we are about to con- 
sider ia designed to remove, as far as poaaible, every arbitrary 
consideration, and to furnish a set of principles which shall 
always guide us to the most probable results. The conscientious 
observer, having taken every care in his observation, will set it 
down, however discrepant it may appear to him, as a portion of 
the testimony collected, out of wliioh the truth, or the nearest 
approximation to it, is to be sifted. 

Admitting, therefore, that the observations give ns the best, 
as indeed the only, information we can obtain respecting the 
desired quantities, we must find a system of corrections which 
shall not only produce the desired accordance, but which shall 
also be the most probable corrections, and further be rendered moat 
probable li^ tliese observations themselves. 

TEE ARITBMEnCAL MEAN. 

4. In order to discover a principle which may serve as a basis 
for the investigation, let us examine first the ease of direct ob- 
servations made for the purpose of determining a single-unknown 
quantity. 

Let the quantity to be determined by direct observation be 
denoted by x. (Suppose, for example, to fix our ideas, that this 
quantity is the linear distance between two fixed terrestrial 
points.) If but one measure of x is taken and the result is a, 
we must accept as the only and, therefore, the most probable 
value, X = a. Let a second observation, takeu under the same 
or precisely equivalent circumstances, and with the same degree 
of care, so that tliere is no reason for supposing it to be more in 
error than the first, give the value b. Then, since there is no 
reason for preferring one obser\'ation to the other, the value of 
X must be so taken that the differences x — a, x — b shall be 
uumerically equal ; and this givee 

x = \ia -\-h) 

This result must be regarded as the only one that can be inferred 
from the two observations cousiateutly with our definition of 
accidental errors; for positive and negative accidental errors of 
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equal absolute magnitude are to be regarded as equal errors and 
as equally probable, since, from the care bestowed on the obser- 
vations and the supposed similarity of the circumstances under 
which they are made, there is no reason a priori for assuming 
either a positive or a negative error to be the more probable. 

Now let a third observation be added, giving the value c. 
Since the three observations are of equal reliability, or, as we 
shall hereafter say, of equal weighty we must so combine a, 6, and 
c that each shall have a like influence upon the result ; in other 
words, X must be a symmetrical function of a, 6, and c, K we 
first consider a and b alone, then a and o, then 6 and o, we shall 
find the values 

with each of which the additional observation c, 6, or a is to be 
combined. Each combination must result in the same sym- 
metrical function, which, whatever it may be, can be denoted by 
the functional symbol o^* We must, therefore, have 

a? = 4' [i (^ + b\ o] 

= 4 [i (^ + c)> «] 
Introducing the sum of a, 6, and c, or putting 

these become 

= 4 [i (« — «)» «] = 4 [«, «] 

But 8 is already a symmetrical function of a, 6, and (?, and there- 
fore these equations cannot all result in the same symmetrical 
function unless c, 6, a, in the respective developments of the 
functions, disappear and leave only s. Hence we must have 

a; = 4(5) 

Now, to determine o^? ^^ observe that, as it must be general, 
its nature may be learned from any special but known case. 
Such a case is that in which the three observations give three 
equal values, or a = 6 = c; and in that case we have, as the 
only value, x = a, or 

a = 4(3 a) 



and, consequently, the symbol ■^ Bigaifies here the div 
Hence, generally, 



In the same manner, if it had been previously shown that for 
m equally good observations the most probable value ia 



it would follow that for an additional observation p we must 
have 

^^ a + b + c+ .... + n+p 
m + 1 

for, putting s = a + b + c + . . . -\- n + p, we shall have 



3; = + |— (s— P).pJ = 



■V [s, J>] = 4 is), &c. 



t we have shown that the form is true for three observed 
values: hence, it is true for four; and since it is true for four 
values it is true for five ; and thus generally for any number.* 

The principle here demonstrated, that the arithmetical mean 
of a number of equally good observations is the most probable 
value of the observed quantity, is that which has been universally 
adopted as the most simple and obvious, and might well be 
received as axiomatic. The above demonstration is chiefly 
valuable as exhibiting somewhat more clearly the nature of the 
assumption that underlies the principle, which is that, under 
strictly similar circumstances, positive and negative errors of the 
aame absolute amount are equally probable. 

5. If now «', n", n'" w'"* aro the m. obser\'ed values of a 

required quantity x, and if x^ denotes their arithmetical mean, 
the assumption of 3:^ as the most probable value of x gives 
n'^x^, n" — x^, n'" — Xt„ &c., as the most probable system of cor- 
rections (subtractive from the observed values) whieh produce 
the required accordance. But the equation 
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may also be put under the form 

that is, the algebraic sum of the corrections is zero. 

This is, however, not the only characteristic of the system of 
corrections resulting from the use of the arithmetical mean. Let 
us examine the sum of the squares of the corrections. For 
brevity, let us denote the corrections, or, as they will be here- 
after called, the residualSy by the symbol v: so that 

i/ = n' — rco, v" = n" — rC(^ t/" = n'" — Xq, &c. 

and also denote the sums of quantities of the same kind by 
enclosing the common symbol in rectangular brackets : so that * 

[r] = 1/ + i/' + t?'" + &c. 

[rv] = t/u' + ly'v" + t;"V" + &c. 

a notation usually employed throughout the method of least 

squares. We have 

M = (2) 

and 

[w] = (n' - x,y + in" - z,)' + (»'" -x^r^ 

== [nn] — 2 [»] a^ + ""^o* 

But since we have also 

[n] 

•X/a — — — — 

m 
this equation becomes 



= [wn] — 






(3) 



m 

Let x^ be any assumed value of x, giving the residuals 

t?j = n' — x^ Vj = n" — iCj V, = n'" — x^, &c. 

then, as above, 

[ Vj = [w**] — 2 [n] ^Tj + mxj* 
Substituting in this the value of [ww] given by (3), we find 

fnV 






(4) 
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This equation determineB the aum of the squares of the residuals 
for any assumed value of x. Since the last term is always posi- 
tive, we see that this sum for any value of x differing from the 
arithmetical mean x^ is always greater than [ru]. Henee it is a 
second characteristic of the arithmetical mean, that it makes the 
sum of the squares of the residuals a minimum. 

6. Obeervations may be not only direct, that is, made directly 
upon the quantity to be determined, but also indirect, that is, made 
upon some quantity which is a function of one or more quanti- 
ties to be determined. Indeed, the greater part of the observa- 
tions ia astronomy, aud in physical science generally, belong to 

the latter class. Thus, let x,y,z be the quantities to be 

determined, and Ma. function of them denoted by/, or 

M = f{x,y,z....) (5) 

and let us suppose an observation to be made upon the value of 
M. We then have but a single equation between x,y, z,. , . and 
the observed quantity M, aud the problem ia as yet indetermi- 
nate. Various syatems of values may be found to satisfy the 
equation, either exactly or approximately. Let us, however, sup- 
pose that the most probable systein (as yet unknown) is expressed 
hy x = p,y = q, z = r...., and let the value of the function, 
when these values are substituted in it, be denoted by F, or put 

V = fip,q,T....) (6) 

then M — V^is the residual error of the observation. In like 
manner, if a number of observations of the same kind be taken, 
in which the observed quantities M', M", M"' . . . are functions 

determined by the same elements p, q, r, , and if V, V", 

V" are the values of these functions when p, q, r. ... are 

substituted in them, then M' — V, M" — V", M'" — V" 

are the residual errors of the observations. If there are /* 
unknown quantities and also n observations, and no more, there 
will be ji equations between the known and unknown quantities, 
which will fully determine the values of these unknown quanti- 
ties: BO tliat the probable values p, q, r....are, in tliat case, 
those determinate values which exactly satisfy all the equations, 
and, consequently, reduce every one of the residuals M' — V, 
M" — V", &c. to zero. But, if there are more than ft observations, 
the determinate values found from fi equations alone will not 
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necessarily satisfy the remaining equations, in consequence of 
accidental errors in the observations. The problem, then, is to 
determine from all t?ie observations^ or from all the equations^ the 
most probable system of values of the unknown quxintitieSj or, which is 
the same thing, the most probable system of residual errors. In the 
case of direct observations, we have seen that the most probable 
value of the unknown quantity was that which made the alge- 
braic sum of the residuals zero ; but this principle followed from 
taking the arithmetical mean of the sam£ quantity, and is ob- 
viously inapplicable in the present case. The second principle, 
that the most probable value is that which makes the sum of the 
squares of the residuals a minimum, is of a more general 
character, and might be assumed at once, as at least a plausible 
principle, to serve as the basis of the solution of our problem ; 
but it will be more satisfactory to justify its adoption by the 
calculus of probabilities. 

THB PROBABILITY CURVE. 

7. Although accidental errors would seem at first sight to be 
of a capricious and irregular nature which would exclude them 
from the domain of mathematics, yet, upon examination from 
theoretical considerations, confirmed, as will be shown, by expe- 
rience, we shall find that they are subject to remarkably precise 
laws. In the first place, we remark that they are subject to the 
following fundamental laws : 1st. Errors in excess and in defect 
— i.e. positive and negative, but of equal absolute value — are 
equally probable, and in a large number of observations are 
equally frequent. 2d. In every species of observations, there is 
a limit of error which the greatest accidental errors do not 
exceed : thus, if I denotes the absolute magnitude of this limit, 
all the positive errors are comprised between and + Z, and all 
the negative errors between and — ?, and, consequently, all the 
errors are distributed over the interval 21. 3d. The errors are 
not distributed uniformly over this interval 2Z, but the smaller 
errors are more frequent than the larger ones. 

Thus the frequency of an error of a given magnitude may be 
regarded as a function of the error itself: so that, if we denote 
an error of a certain magnitude by J, and its relative frequency 
in a given large number of observations by ^J, this function 
should obtain its maximum value for J = 0, and become zero 
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when J = ± ^. If, then, 
J by y, or put 



re denote the probability* of an error 

y = V^ (7) 

we may regard this as the equation of a curve, taking J as the 
abeciesa and y as the ordinate. The nature of this curve will be 
accurately defined when we have discovered the form of the 
fuuetiou fJ, but we can see iu advance that a curve such as 
Fig. A is required to satisfy the couditioua already imposed upon 




! function. For its maximum ordinate must correspond to 
J = ; it must be symmetrical with reference to the axis of y, 
since equal errors with opposite signs have equal probabilities; 
and it must approach very near to tlie axis of ahscisase for values 
of J near the extreme limits, although the impoasihility of aa- 
eigning such extreme limits of error with precision must prevent 
us from fixing the poiut at which the curve will finally meet the 



8. The number of possible errors in any class of obsei-vations 
is, strictly speaking, finite ; for there ia always a limit of accuracy 
to the observations, even when we employ the moat refined 
instruments, in consequence of which there is a numerical suc- 
cession in our results. Thus, if 1" is the smallest measure in a 
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given case, the possible errors, arranged in their order of magni- 
tude, can only differ by 1" or an integral number of seconds. 
Hence, our geometrical representation should strictly consist of 
a number of isolated points ; but, as these points will be more 
and more nearly represented by a continuous curve as we increase 
the accuracy of the observations, and thus diminish the intervals 
between the successive ordinates, we may, without hesitation, 
adopt such a continuous curve as expressing the law of error. 
We shall, therefore, regard J as a continuous variable, and <pA 
as a continuous function of it. 

Now, by the theory of probabilities, if ^ J, ^ J', ^ J" 

are the respective probabilities of all the possible errors J, J', 
J" we have* 

when the number of possible errors is finite. But the assumed 
continuity of our curve require? that we consider the difference 
between successive values of J as infinitesimal, and thus the 
number of values of ipA is infinite, and the probability of any 
one of these errors is an infinitesimal. To meet this difficulty, 
let us observe that if a finite series of errors J, J', J" .... be ex- 
pressed in the smallest unit employed in the observations, these 
errors, arranged in the order of their magnitude, will be a series 
of consecutive integral numbers ; the probability of the error J 
may be regarded as the same as the probability that the error 
falls between J and J + 1 ; and the probability of an error be- 
tween J and J + I will be the sum of the probabilities of the 
errors J, J + 1> ^ + 2, ^ + {i — 1). If Hs small, the pro- 
bability of each of the errors from J to J + f will be nearly the 
same as that of J : so that their sum will differ but little from 
i<pj. As the interval between the successive errors diminishes, 
this expression becomes more accurate ; and hence when we take 
rfJ, the infinitesimal, instead of i, we have ipA . dd as the rigorous 
expression of the probability that an error falls between J and 
J + dJ. Hence, it follows, in general, that the probability that 
an error falls between any given limits a and 6 is the sum of all 

* For if there are n errors equal to J, n' equal to J', &€., and the whole number 

of errors is m, the probabilities of the errors are respectively 0J = — , 0J'=: — , &c., 

n 4- n' 4- » . , . m ^ ^ 

and the sum of these is — -J- — ^ * ' ' ' = — = 1. 

m tn 
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the elements of the form <pJ . dJ between these limits, or the 
Integi 



ntegral, taken between the extreme limits of error, and 
thus embracing all the possible errors, will be 



/:, 



,j.dj^\ 



We liave heretofore assumed that the function (pA is to bo zei-o 
for A = ± I. It must also be added that, since the probability 
of any error greater than ± Ms also zei-o, we should have to 
determine this function in such a manner that it would be zero 
for all values of J from + Ito + x and from — Ito — x. The 
obvious impossibility of determining snch a function leads ua 
to extend the limits ± I to ± <x, aud to take 



I to ex ti 



/_: 



f.J.rfJ=l (8) 



This will evidently be allowable if the integral taken from 
±lto±:'X! is so small as to be practically insignificant. Besides, 
the extreme limits of error can never be fixed with precision, aud 
it will suffice if the function ipJ is such that it becomesvery small 
for tliose errors which are regarded as very large, 

9. Returning now to the general case of indirect observations, 
Art. 6, in which we suppose a quantity M = f(x, i/, z, . . . .) to he 

observed, let J, J', J" be the errors of the several observed 

values of M, and yj J, p J', p J" their respective probabilities ; 

then, the probability that these errors occur at the same time in 
the given series being denoted by P, we have, by a theorem of 
the calculus of probabilities,* 

P= f.J.pJ'.y.J" (9) 

The moat probable system of values of the unkiio^vn quantities 

* If a Binglo Bction of a cause can produce [ha effcoO u, n', a" wilh Ihe ro- 

Bpeeli»e probabiUriee p, p', p", Ihe prpbnlJilily that Iwo Bueceasive iodependeol 

actions of the cnaae will produce the effecia o sni a' h pp': aiiJ aimilsrlj for any 
number of etfeots. Thu!. if an urn contains 2 while balls, 3 red ones, and 6 blnck 
one*, the probability thai in Iwo auoceeaive drawing! (Ihe original number of bftUs 
being (he same al each drawing] one ball will be white and the other red >biV X i'^* 
„ II,— 11 
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x, y, ^ . . . . will be that which makes the probability P a maxi- 
mum. Consequently, since x, y, ^ . . . . are here supposed to be 
independent,* the derivative of P relatively to each of these 
variables must be equal to zero ; or, since log P varies with P, 
the derivatives of log P must satisfy this condition, and we shall 

have 

1 dP ^ \ dP ^ . 

-..--=0, —.—- = 0, &c. 

P dx P dy ' 

which, since 



log P = log f J + log f J' + log f J" + ...., 

give the equations 

dx dx dx 

^'■d^ + ^^-dH + ^'^W^ > (10) 

az dz az 

kc. &c. 

in which we have put 

fj.dd 



^^=^rT-. (11) 



The number of equations in (10) being the same as that of the 
unknown quantities, these equations will serve to determine the 
unknown quantities when we have discovered the value of the 
function ^' J, as will be shown hereafter. 

Since the functions ipA and ^' J are supposed to be general, and 
therefore applicable whatever the number of unknown quanti- 
ties, we may determine them by an examination of the special 
case in which there is but one unknown quantity, or that in 
which the observed values J!f, J!f' , J!/". . . . belong to the same 
quantity. In that case, the hypothesis that x is the value of this 
quantity gives the errors 

A = M — x, A'=M' — x, A" = M" — x 

* That is, subject to no restrictions except that they shaU satisfy the obserrations, 
or the equations if =/(z, y, s, . . . .). For the case of ** conditioned" obseryationB, 
see Art. 58 of this Appendix. 
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whence 

dJ dA' dJ'' , 

dx dx dx 
and the first equation of (10) becomes 

/ (Jf — x) + ^ (JT — ^) + sp' (^" — ^) + = (12) 

This being general for any number m of observations, and for 

any observed values J!f, M\ iHf ". . . ., let us suppose the special 

case 

M' = 2r* =M--mN 

Since the arithmetical mean of the observed quantities is here 
the most probable value of x, we have 

x = l(jf+jf'+jr'+ ) 

= l[Jf+(m-l)(Jf-m2r)] 



whence 



= Jf — (m — 1) iT 

Jf— x = (m — l)if 

M'-^x = M''-^x == — JV^ 



and, consequently, (12) becomes 

SP' [(m — 1) JVT] 4. (m — l)/(— JVr)=rO 

y^[(m~l)iy] ^ y^(-JO 

That is, for all values of m, and therefore for all values of (m — l)iV, 

we have — 7^^ ^^ ^^ equal to the same quantity ^— ^^ — ri~- 

' ©'J 

Hence we have generally — equal to a constant quantity, and, 

denoting this constant by A, we have 
or, by (11), 

Integrating, 

log f J == } A: J* + log x 
whence 

in which 6 is the base of the Napierian system of logarithms. 
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Since ipA must decrease as J increases, \k must be essentially 
negative : representing it, therefore, by — A*, our function becomes 



— V fi^hhAk^ 



To determine the constant x, let this value be substituted in (8), 
which gives 

Putting ^- 

t = AJ (13) 

this integral becomes 

The known value of the definite integral in the first member is 
|/;r (see Vol. I. p. 153) ; whence 



X = 



V^ 



and the complete expression of <pA becomes 






(14) 



The constant h must depend upon the nature of the observa- 
tions, and will be particularly examined hereafter. If we here 
take it as the unit of abscissae in the curve of probability, the 
equation (7) becomes 

1 



>— AA 



by which the curve may be constructed. The values of y for a 
few values of J are as follows: 



A 


y 


0.0 


0.5642 


0.2 


0.5421 


0.4 


0.4808 


0.6 


0.3936 


0.8 


0,2975 


1.0 


0.2076 


1.2 


0.1837 


1.4 


0.0795 


1.6 


0.0486 



Diff. 



.0221 
.0613 
.0872 
.0961 
.0899 
.0739 
.0542 
.0359 





J 


y 


1.6 


0.0486 




1.8 


0.0221 




2.0 


0.0103 




2.2 


0.0045 




2.4 


0.0018 




2.6 


0.0007 




2.8 


0.0002 




3.0 


0.0001 




00 


0.0000 



Biff. 



.0215 
.0118 
.0058 
.0027 
.0011 
.0005 
.0001 
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The curve, Fig. A, in Art. 7, ie constnicted from this table ; but, 
to exhibit its character more dietiactly, the scale of the ordiiiatea 
is four times that of the absciasfe (which, indeed, corresponds to 
the case oi k = 2). "We see that the curve approaches very near 
to the axia for moderate values of J, and that the assumption of 
± cc instead of finite limits of J can involve no practical error. 
It is evident that the axis J'X is an asymptote to the curve. 

The differences in the above table indicate that the curve 
approaches the axis most rapidly at a point whose abscissa is 
between 0.6 and 0.8, The exact position of tliis point, which 
is a point of inflexion, is found by putting the second differen- 
tial coefficient of j/ equal to zero, which gives 

£, r— 44 J ff — 44 n 



-V2- ■■ 

The ordinate Mj/i is drawn at this point. We shall have o 
sion to refer to it again hereafter. 



THE MEASURE OF PRECISION. 

10. The constant h requires special consideration. Since the 
exponent of e in (14) must be an abstract nnmber, r must be a 
concrete quantity of the same kind as J. In a class of observa- 
tions in which J is small for a given probability fJ, j will be 
small, and h will be large. Thus, k will be the greater the more 
precise the nature of the observations, and is, therefore, called by 
Gauss the measure of precision. If in one syetem of observa- 
tions the probability of an error J is expressed by 



and in another, more or less precise, by 



li and in 



V" 



tlie probability that in one observation of the first system the 
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error committed will be comprieed between the Umita - d and 
+ i will be expressed by the integral 



X 






and, in like manner, the probability that the error of an observa- 
tion in the second system will be comprised between - d' and 
+ i' will be expressed by 

X+y A' 

These integrals are evidently equal when we have h8 = A'^'. If, 
for example, we have A' = 2A, the integrals will be equal when 
3 == 2^'; that is, the double error will be committed in the first 
system with the same probability as the simple error in the 
second, or, in the usual mode of expression, the second system 
will be twice as precise as the first. "We shall presently see how 
the value of h can be found for any given observations* 

THE METHOD OF LEAST SQUARES. 

11. The preceding discussion leads directly to important prac- 
tical results. We have seen (Art. 9) that to find the most probable 

values of a:,y, ^ . . . . from the observed values of J!f =/(a:,y, ^, ) 

we are to render the probabilityP= yJ.^J'. yJ". ... a maxi- 
mum, that is, by (14), 

must be a maximum; and this requires that the quantity 
J J + J' J' + J" J" + • . . • should be a minimum. Thus, the prin- 
ciple that the most probable values of (he unknown quantities are those 
which make the sum of the squares of the residual errors a minimum^ is 
not limited to the case of direct observations, but is entirely 
general. 

The principle is readily extended to observations of unequal 
precision. For if the degree of precision of the observations 
M^ M'^ J!f ". ... be respectively A, A', A". . . ., and we compare 
these observed quantities with the values F", F', T". . . ., computed 
with the most probable values of a:,y,^. . . ., whereby we obtain 
the residual errors M — V^= J, M* — y == J'. . . ., it is the same 
thing as if we had taken observations of equal precision (repre- 
sented by 1) upon the quantities Ailf, A'ilf' , A"JK'". . . ., and had 
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compared them with the computed quantities h V, h' V, h" V" . 

whereby we should have found the errors hM — kV ^^ hJ, 

h'M' — h.'V'^= h'A' , iu which case we should have to reduce 

to a minimum the quantity 

A' J' + ft" J" + A'" J'" + . . . . 

that is, each error beivg wuhiplied by its measure of precision, nml 
therehy reduced to the same degree of 'precision, the. sum of the squares 
of the reduced errors must be a minimum. 

In what precedes is involved the whole theory of the method 
of least squares. I proceed to develop its practical features. 

THE PROBABLE ERROE, 

12. From the preceding articles it follows that the prohability 
that the error of au observation falls between J and J -\- dJ \s 
expressed by 

and the probability that it falls behveen tlie limits and a is 
expressed by 

and this integral expresses the number of errors that we should 
expect to find between the limits and a when the whole num- 
ber of errors is put — 1 [equation (8)]. If we put / — hA, the 
integral takes the form 






The whole number of errors, both positive and negative, whose 
UQmerical magnitude falls between the given limits is twice this 
integral, or 

The value of this integral (which may be computed by the 
methods of Vol. L Art. 113) is given in Table IX. The number 
of errors between any two given limits will be found by taking 
the difference between the tabular numbers corresponding to 
these limits. Since the total number of erroi-s is taken as unity 
in the table, the required number of errors in any particular case 
is to be found by multiplying the tabular numbers by the actual 
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number of observations. Thus, if there are 1000 observations, 
we find that 



between 

u 
u 
l( 
ti 



= and t = 0.5 there are 520 errors. 
= 0.5 " t = 1.0 « " 322 " 
= 1.0 " t = 1.5 " *' 123 " 
= 1.5 " t = 2.0 " ^* 29 " 
= 2.0 " <= 00 « « 5 " 



13. The degrees of precision of different series of observations 
may be compared together either by comparing the values of A, 
or by comparing the errors which are committed with equal 
facility in the two systems. The errors to be compared must 
occupy in the two systems a like position in relation to the ex- 
treme errors, and we may select for this purpose in each system 
the error which occupies the middle place in the series of errors arranged 
in the order of their magnitude^ so that the number of errors which are 
less than this assumed error is the same as the number of errors which 
exceed it. The error which satisfies this condition is that for 
which the value of the integral (16) is 0.5. Denoting the cor- 
responding value of t by />, we find, by interpolation from Table 

IX., 

P = 0.47694 
and we have 

2 /•p 1 

If then we denote by r the error which, in any system of obser- 
vations whose degree of precision is A, corresponds to the value 
t = py or put 

p = hr h = ^ (18) 

there will be a probability of J that the error of any single obser- 
vation in that system will be less than r, and the same proba- 
bility that it will be greater than r ; which is sometimes expressed 
by saying that it is an even wager that the error trnU be less than r. 
Hence r is called the probable error. 

We may, therefore, compare different series of observations 
by comparing their probable errors, their degrees of precision 
being, by (18), inversely proportional to these errors. 

14. In order to apply Table IX. in determining the number 
of errors in a given class of observations, we must know the 
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measure of precision A, or the probable error r.- thus, if we 
wieh the number of errors less than a, we enter the table with 
the argument t = ah, or ( ^ — 

For greater convenience, we cam employ Table IX.A, which 
^ves the same function with the argument -. For example, if 
there are 1000 observations whose probable error ia ?■ ^ 2", 
and we wish to know the number of errors less than a = 1", we 
take from Table EX.A, with the argument - = 0,5, tlie number 
0.26407, which multiplied by lOOO gives 264 aa the req^uired 
number. 

The following example from the Fundamenta Astranomi<B of 
Bebsel will serve to show how far the preceding theory is sus- 
tained by experience. In 470 observations made by Biladlby 
upon the right ascension of Sinus and AUair, Bessel foujid the 
probable error of a single observation to be 

r = 0".2637 

Hence, for the number of errors leas than 0".l the argument of 
Table IX.A will be ^|i- = 0.3792; and for 0."2, 0".3, &c., the 
3 multiples of 0.3792. Thus, we find from the table 



for 0".l with arg 


0.8792 the 


nomber 0.20187 


" .2 


0.7584 


0.39102 


" .3 


1.1376 


0..55710 


" A 


1.5168 


" 0.69372 


" .6 


1.8960 


« 0.79904 


" .6 


2.2752 


" 0.87511 


■' .7 


2.6544 


' 0.92661 


" .8 " 


3.0336 


" 0.95926 


" .9 


3,4128 


" 0.97866 


" 1 .0 


3.7920 


" 0.98946 
" 1.00000 



Subtracting each number from the following one, and muttiply- 
ing the remainder by 470, the number of observations, there were 
found 
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Between 


No. of errors 
by the theory. 


No. of errors by 
ezperienoe. 


0".0 and 0".l 


96 


94 


.1 " .2 


89 


88 


.2 " .8 


78 


78 


.8 " .4 


64 


58 


.4 " .6 


50 


51 


.5 " .6 


86 


86 


.6 " .7 


24 


26 


.7 " .8 


15 


14 


.8 " .9 


9 


10 


.9 " 1 .0 


5 


7 


over 1 .0 


5 


8 



The agreement between the theory and experience, though 
not absolute, is remarkably close. The number of large errors 
by experience exceeds that given by the theory, and this has 
been found in other cases of a similar kind; which shows at least 
that the extension of the limits of error to dz oo has not intro- 
duced any error. The discrepancy rather indicates a source of 
error of an abnormal character, and calls for some criterion by 
which such abnormal observations may be excluded from our 
discussions and not permitted to vitiate our results. Such a 
criterion has been proposed by Prof. Peirce, and will be con- 
sidered hereafter. 



THE MEAN OF THE ERRORS, AND THE MEAN ERROR. 

15. The selection of the probable error as the term of com- 
parison between different series of observations is arbitrary, 
although it seems to be naturally designated by its middle posi- 
tion in the series of errors. There are two other errors which 
have been used for the same purpose. 

The first is the m^an of the errors^ these being all taken with 
the positive sign. In order to find its relation to the probable 
error, let us first consider a finite series of errors 

with the respective probabilities 



2a 



2a' 



2a" 
m ' 



1 



METHOD OF LEAST SQUARES. 491 



I SO that in m observations there will be 2a errors (tnimerically) 
equal to J, 2a' equal to J', &c., the probability of a positive error 
J being — . The mean of all these errors, each being repeated 
a number of timea proportional to its probability, is 

When the number of errors is infinite, the probability of an 
error J is to be understood as the probability that it falls 
between J and J -f dJ, which is ^J.dJ (Art. 8), and the above 
formula for the mean of the errors becomea the sum of an infi- 
nite number of terms of the form 2JfJ.dJ. Hence, putting 

i; =; the mean of the errors, 



=/; 



or, by (18), 



" fev-ir 



C19) 



^^^1 



^ 0.8453 1) 



Another error, very commonly employed in expressing the 
precision of observations, ia that which has received the appella- 
tion of the mean error {der mittlere Fehler of the Germans), which 
is not to be confounded with the above mean of the errors. Its 
definition is, the error Ike square of which, is the mean of (He squares of 
all the errors. Hence, putting 



(21) 





e ^= the mean error 


we have 






■■=r."7j^'«-""" 


or, by (18), 






t = -^ = 1.4826 r 

r^ 0.6745 c 



(22) 



When we put ft — 1, we have s — vi- The mean error is, 
therefore, tlie abBcissa of the point of inflection of the curve of 
probability (Art. 9). In the figure, p. 479, OM is the mean error, 
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OP the probable error, OE the mean of the errors, and -Mm, i^, 
JEcy their respective probabilities. 



THB PROBABLE ERROR OF THE ABITHMBTICAL MEAN. 

16. The error above denoted by r is the probable error of any 
one of the observed values of the unknown quantity x. We are 
next to determine the relation between this and the probable 
error r^ of the arithmetical mean of these values. 

If J, J', J" . . . . are the errors of the observed values, the 
most probable value of x is that which renders the probability 

P = A«»7r— i«*^— AA(AA + A'A' + A"A''+...) 

a maximum (Art. 11), and, consequently, the sum JJ + J' J' 
+ .... a minimum. But this sum is rendered a minimum by 
the assumption of the arithmetical mean x^ as the most probable 
value (Art. 5), and hence the quantity P expresses the probability 
of the arithmetical mean if J, J', J" . . . . are the errors of the 
observations when compared with this mean. The probability 
of any other value of x, as ar^ + ^, will be 

Since [J] = J + J' + J" + ....== (Art. 5), and [JJ] = mu 
(Art. 15), this expression may be put imder the form 

and at the same time we have 

so that 

jP: jP' = 1 le—vthhii 

that is, the probability of the error zero in the arithmetical mean 
is to that of the error <J as 1 : e'-^'*^. For a single observation, 
the probability of the error zero is to that of the error 8 as 
l:e— **". Hence the measure of precision (Art 10) of the 
single observation being A, that of the arithmetical mean of m 
such observations is hi/m; from which follows the importwnt 



I 



^/ f ^ 



'MLfU/r^'e 
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theorem that ike ■precision of the mean of a number of obsercalions 
increases as the square root of their number,* 

If, then, r is the probable error of a single obaer\'atiou, aud r 
that of the arithmetical mean, we must have 



l/m 



(23) 



au(i from the conBtant relation between the mean and the proba- 
ble error (22), ^ 

^. = -~ C24) 



DETERMINATION OF THE MEAN AND PROBABLE ERRORS OP GIVEN 
OBSERVATIONS. 

17. The principles now explained will enable U8 to determine 
the mean errors of any given series of directly observed quanti- 
ties. Let n, n', n" be the observed values; x, tiieir arith- 
metical mean ; v, v', c" . . . . the residuals found by subtracting 
Xn from each observed value : so that 

V = n — x^, xf =^n' — Xa, "" =^ n" — x^, &c. 

If Xa were certainly the true value of x, so that v, v', v" were 

the actual or (as we may say) the true errors, and, consequently, 

identical with J, J', J" , we should have, according to the 

l^jhore, mee = [JJ] — [ctj], and hence 

and this must always ^ve a close approximation to the value of e. 

But the relation mse = [JJ] was deduced from a consideration 

of an infinite aeries of errors which would reduce the mean 

error of Xg to an infinitesimal, according to the principles assumed, 

; and thus make v, v', v" .. , . identical with J, J', J" . . . A better 

i approximation to the value of t, where the series is limited, is to 

[ be obtained by considering the mean error of x^ itself, and conse- 

J quently, also, the mean errors of the residuals ;;, v', v" If 

[, then we suppose the true value of x to be x^-f S, we shall have 
the true eri-ors 



■-m 
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whence, observing that [r] = 0, 

= [^ + ^^' 

Thus the approximate value nui = [vv] requires the correction 
m8*, the value of which depends upon the value we may ascribe 
to d. As the best approximation, we may assume it to be the 
mean error Cq : so that, by (24), 

• m 

which gives 

nut = [yv] + «« 
whence 



n 



=S -4B^) <»> 



and consequently, also, by (22), 



r=*q 



\/(S) '=«■"« w 



Thus from the actual residuals the mean and the probable error 
_pf a single observed value are found. Hence, by (28) and (24), 
the mean and probable errors of the arithmeti^l mean will be 
found by the formulae 






«0 



=V(i|^) '.=W(i^) w 



Example. — ^Let us take the following measures of the outer 

diameter of Saturn's ring observed by Bessbl at the Eonigsbei^ 

Qbservatory with the heliometer, in the years 1829-1831.* The 

measures, denoted by n, are all reduced to the mean distance of 

* Saturn from the sun, and are here assumed to have the same 

degree of precision. 

Il l . . . . 

* Attrw. Naeh., Vol. XII. p. 169. 
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38".91 


— 0".40 


0.1600 




39".41 


+ 0' .10 


0.0100 




39 .82 


+ .01 


.0001 




39 .40 


+ .09 


.0081 




38 .93 


-0 .38 


.1444 




39 .36 


+ .06 


.0025 




39 .31 


.00 


.0000 




39 .20 


-0 .11 


.0121 




39 .17 


-0 .14 


.0196 




39 .42 


+ .11 


.0121 




39 .04 


-0 .27 


.0729 




39 .30 


-0 .01 


.0001 




39 .67 


+ .26 


.0676 




39 .41 


+ .10 


.0100 




39 .46 


+ .15 


.0225 




39 .43 


+ .12 


.0141 




39 .30 


-0 .01 


.0001 




39 .43 


+ .12 


.0144 




39 .03 


-0 .28 


.0784 




39 .36 


+ .06 


.0025 




39 .35 


+ .04 


.0016 




39 .02 


-0 .29 


.0841 




39 .25 


-0 .06 


.0036 




39 .01 


-0 .30 


.0900 




39 .14 


-0 .17 


.0289 




38 .86 


-0 .45 


.2025 




39 .47 


+ .16 


.0266 




39 .61 


+ .20 


.0400 




39 .29 


-0 .02 


.0004 




39 .21 


^0 .10 


.0100 




39 .32 


+ .01 


.0001 




39 .17 


-0 .14 


.0196 




39 .40 


+ .09 


.0081 




39 .60 


+ .29 


.0841 




39 .33 


+ .02 


.0004 




39 .54 


+ .23 


.0629 




39 .28 


— .03 


.0009 




39 .45 


+ .14 


.0196 




89 .62 


+ .31 


.0961 




39 .72 


+ .41 


.1681 




^0 - 


= 39 .308 


M = 


= 1.5884 




Hence, eince m = 40, we have, by (25) and (26), 


^ 'W(f^)-°"-^ a 


^^P r 0".302 X 0.6745 = 0".136 ^^H 


and couBequently, by (23) and (24), or (27), ^^^H 


0" 202 0" 136 ^^H 


"■■'"■' 0".032, r, — 0".022 


" 1/(40) ' • vm 


That is, the probable error of a eingle obaervationjvae 0",136, 


aud that of the final reault x^ = 39".308 was only 0".022. 


18. The preceding method of finding the probable error from 


the squares of the residuals is that which is moat commonly 


employed ; but when the number of observations is very great, 


' it is desirable to abridge the labor, if possible. A eufflcient 


approximation can be obtained by the use of the first powers of 


the residuals as follows. 


The number of observations being very great, we ahall pro- 


bably have as many positive as negative residuals. If r', v", 


L _y 


l^^^^^^^^l^^^^ 
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t?'" . . . are the positive and Vj, r,, r, . . . the negative residuals, 
and if the true value of x is x^-i- 8, the true errors will be 

v' — d, r" — d, v'" —d , and — t?i — *, — i?, — *, — v, — ^, 

If they are all taken loith the positive sign only^ the errors are, there- 
fore, 

1/ — dj v" — d, I?"' — ^, . . . . and Vj + ^, t>, + ^, v, + ^, 



the mean of which, upon the hypothesis of an equal number of 
positive and negative residuals, is the same as that of the series 

Hence, denoting the sum of the numerical values of the residuals 
by [r], and the mean of the actual errors by jy, as in Art. 15, we 
have 



and hence, by (20), 






r = 0.8453 ^ (28) 



and consequently, also, by (22), 

c = 1.2533 ^ (29) 

In the example of the preceding article we find the mean of the 
residuals taken with the positive sign to be 0".1565, which by 
(28) gives r = 0'M665 X 0.8453 = 0".131, which is perhaps a 
sufficient approximation to the value found above. In this 
example, however, we have 22 positive residuals, 17 negative 
ones, and 1 zero: so that the hypothesis upon which the formula 
(28) was founded is not strictly applicable. In a larger number 
of observations we should expect a closer agreement with the 
hypothesis, and more accordant results. 

We may, however, employ the first powers of the residuals 
more strictly according to the theory of probabilities. In a 
limited series each residual is to be regarded as liable to a pro- 
bable error r', and their mean is to be regarded as the mean of 
the errors of the residuals themselves, rather than as the mean 
of the errors of the observations. Hence the formula 



r' = 0.8453 

m 



\ 
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gives the probable error of a residual. The relation between 
r' and r (= the probable error of an observed quantity 7i) may be 
found as follows. Each observed n may be supposed to be the 
result of observing the mean quantity Xq increased by an ob- 
served error v. The probable error of n = x^ + v is, therefore 
(by a principle hereafter to be proved), 



whence 



r = i/(ro' + = ^/(^' + r'^) 



= ''^^J 



or 



r = 0.8453 ^ (30) 

v/[m(m — 1)] ^ ^ 

which agrees with the formula given by C. A. F. Peters.* Ac- 
cording to this formula, we find in the above example r = 0".133. 

DETERMINATION OF THE MEAN AND PROBABLE ERRORS OF FUNCTIONS 

OF INDEPENDENT OBSERVED QUANTITIES. 

19. Suppose, first, the most simple function of two independ- 
ent observed quantities x and x^y namely, their sum or difference 

X =^ X ±z x^ 

and let the given mean errors of x and x^ be e and e,. Although 

the number of observations by which x and Xj have been found 

may not be given, we may assume it to have been any large 

number m, and the same for each of the quantities ; the degrees 

of precision of the two series being inversely proportional to e 

and ep The true errors of the assumed observations may be 

assumed to be — 

for Xy J, J', J" 

for X,, J„ J/, J," ...... 

and the errors of X, consequently, 

Denoting the mean error of JTby E^ we have, by the definition, 

mE^ = (J di J,)« + (J' =t J,')' + (-^" ± ^i")' + 

*A8tron, Nach.y Vol. XLIV. p. 82. 
Vol. II.— 32 
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In a great number of observations there must be as many posi- 
tive as negative products of the form J Jj, and such that we shall 
probably have \_JJ^ = ; and since we also have we* = [^4]> 
me^ = [^i-^J, this equation gives 

^« = t« + ei« (31) 

If we have 

JT =^ X :h Xi ±z Xg 

and the mean errors of x, x^, x^ are e, e^ Cj, we have by the pre- 
ceding equation the mean error oi x ± x^ = \/{f + ^i)', and by 
a second application of the same equation, considering x db ajj as 
a single quantity, the mean error of JT will be found by the 
formula 

j&« = e» + ci« + e,« (31*) 

and the same principle may be thus extended to the algebraic 
sum of any number of observed quantities. 
In consequence of the constant relation (22), if r, r^, r, . . . . 

are the probable errors of x, x^, x^ and ^ the probable error 

of ^= X ±1 x^ ± x^ , we shall have 

72* = r« + ri» + r,« + . . . . (32) 

Example 1. — The zenith distance of a star observed in the 
meridian is 

C = 21*^ 17' 20".3 with the mean error e = 2".3 

and the declination of the star is given 

d = 19*^ 30' 14".8 with the mean error ei = 0".8 

Required the mean error JE^ of the latitude of the place of obser- 
vation, found by the formula f = C + ^« We have, by (31), 

E = i/[(2.3)« + (0.8)'] = 2".44 
Hence 

^ = 40*^ 47' 35'M with the mean error E = 2".44 

Example 2. — The latitude of a place has been found with the 
mean error e = 0".25, and the meridian zenith distance of stars 
observed at that place with a certain instrument has been found 
to be subject to the mean error e^ = 0".62 : what is the mean 
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error E of the declinations of the stars deduced by the formula 
5 = y — C? We have 

E = i/[(0.25)« + (0.62)2] ^ o".67" 

20. Let us next consider the function 

X=ax 

and suppose x has been observed with the mean error e, and a is 
a given constant. Every observation of x with the error d= J 
gives JT with the error ±aJ: so that the mean error of Xmust be 

E = 06 

In general, by combining this with the preceding principle, if 

we have 

^ z= ax -{- aiOTi + ^2^2 + . . . . 

and if the mean errors of x, x^, x^ . . . . are e, e^, e^y ? ^^^ -^ 

that of JT, we shall have 

B^ = aV + a^Uy^ + a,U^ + -/^ • = [^'^^ (38) 

and the same form may be used for probable errors. 

Example. — As an example illustrating the application of both 
the preceding principles, suppose that in order to find the rate 
of a chronometer we find at the time t its correction + 12"* 13*.2 
with the mean error 0*.3, and at the time V the correction 
+ 12^ 21*.4 with the same mean error O'.S, and the interval V — t 
= 10 days. The rate in the whole interval is 

12- 21'.4 — 12- 13'.2 = + 8*.2 

with the mean error, according to Art. 19, 

|/[(0.8)« + (0.8)*] = 0'.42 

The mean daily rate is then 

+ ^ = + 0'.82 

with the mean error, according to Art. 20, 

0*42 

_1 = 0'.042 

10 
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. 21. If X, x^y Xj . . . . are the several observed values of the 
same quantity, their arithmetical mean being 

Xq =^ ^ yX -f- Xi — p Xf —(" • • • ) 

tn 

andjf r js..tbe probable error of each 9 b*^^n^Ptir>Pj what is the 
probable error r^ of Xq? By Art. 19, the probable error of the 

sum X + x^ + X2 + . . . . iQ 

l/(r' -f r' + r* + . . .) = ^/(mr^) = ri/m 
and the probable error of — th of the sum is, by Art. 20, 

r, = - X ri/m = —- 
m '^ l/m 

as has been otherwise proved in Art. 16. 

22. Let us now take the general case in which X is any func- 
tion whatever of the observed quantities x, x^, x^, . . . . expressed 
by 

^ J \Xy Xlf X^j » • » *J 

Let the variables be expressed in the form 

UC ■ ' CL "T" 3u y **^\ -~~ vtj *T~" X\ « iCTj ■ dj "l" **'2 I • • • • 

a, ttp ^2 . . . being arbitrarily assumed very nearly equal to 

x, Xi, Xj . . . respectively, and such that x', x/, x^^ may 

be so small that their squares will be insensible. The given 
mean errors e, e^ e^ • • • ^^7 then be regarded as the mean 
errors of x', x/, Xj' . . . . The function X developed by Taylor's 
theorem is 

X=/(a, a„a, ....) + -^^ + ^^1 + ^-^2 + •••• 

and the mean error of JT will be that of the quantity 

dX , r/X , ^JT , 

(/Jl- dx^ dx^ 

. or, by (33), 
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or, if r, Ti, r^ . . . are the probable errors of a:, x^y a:^ . . ., and H 
that of -Z, 

This formula is, indeed, but approximative, since we have 
neglected the terms involving the higher powers in the develop- 
ment of JT; but the mean errors of these small terms will be in- 
sensible if we suppose that the errors e, e^, £3 . . . are so small 
that the differences between the observed values x, x^, x^. .. 
and the true values are of the s am^ ^ gyder as the quantities 
x', x/, x^ .-. ., which will always be the case where proper care 
has been taken to reduce the accidental errors of observation to 
their smallest amount. If the given function is implicit, as 

0=/(X, a:, a?j, ar, ...) 

we should still by differentiation obtain the differential coeffi- 
cients, and then find the mean errof of X by (34). 

Example.— The local apparent time at a place in latitude 
ip = 38° 58' 53" was found (Vol. I. Art. 145) from the sun's 
zenith distance f = 73° 12' 25", when the declination was 
a = — 22° 60' 27", to be < = 2* 47"* 39'.4. What is the probable 
error of this result, supposing the probable errors of the data 
to be — 

Probable error of ^ = r = 0".5 
" " <J = fj = .6 

" *' C = r, = 3 .5 

The formula 

= — cos C + sin v> sin ^ + cos ip cos b cos t 

expresses ^ as an implicit function of f>, d, and f. We find 
(Vol. I. Art. 35) 

dt 1 



d^ 
dt 


cos f> tan A 
1 


dd 
dt 


cos d tan q 
1 


d: 


cos f sin A 



J 
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where A is the azimuth and j the parallactic angle. We find 
from the data A = + 40° 1', q = 32° 51', whence 

^ = _ 1.532, — = 1.680, — = +^.001 
dfp dd dZ 

and the probable error of i is, by (34*) 
R = |/[(0.5 X 1.532)« + (0.6 X 1.680)« + (3.5 X 2.001)*] = r.l2 

or, in seconds of time, 

B = 0-.47 

23. To complete this branch of our subject, it is to be observed 
that the preceding demonstrations apply only to the case wher.e 
the quantities entering intdpiimbination are independent ; but 
when they are merely differdirf; functions of the same observed 
quantities, the above formulae are incomplete. Let us suppose 
that we have Xand X', diflSjyent functions of the same observed 
quantities x, Xp x^, ... ., 



amere 



JT =f (x,x^,x^, ) 

JT' = f\x, x^, x^, ) 

the mean errors of x, x^^ x^ .. . being e, e^, gj • • • • 5 ^^^ ^^^t we 
wish to find the mean error JS of the function. 

If any single observation of x, x^j x^ is aflfected by an error 

^) ^v ^2> • • • • respectively, the corresponding errors in JTand X' 
will be— 

Error in X, J = a^ + ^i^i + ^j^j + 

" ^, J' = a'd + a^d^ + a,M, + . . . . 

in which a, a^, ag . . . . are the differential coefficients of JT, and 

a', a/, Og' the differential coefficients of JT', with reference 

to x^ XjL, ajj, .... The corresponding error in Y will be 

J" =AA + A' A' 

in which A and A' are the differential coefficients of Y with re- 
ference to X and X'. The square of the mean error E will be 
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tlie mean of the squares of all the values of J" which result 
from all the possible values of 5, ^^ dg* • • • 
Substituting the values of J and J', we have 

J" == {Aa + A' a!) d + {Aa^ + A'a^) ^^ + . . . . 

which we may briefly express as follows: -^ 

If the number of values of J" is denoted by m, the mean of all 
the values of J"^ will be 

' m 

In consequence of the various si^|^)f 88^, dd^y &c., the mean 
value of each of these quantities^^Bbe zero ; and the mean 
values of ^, 8^^, &c. are e^, e^*, &c. iRice the formula becomes 
simply 

jE7« = (^Aa + A' ay e» + {Aa^ + ^'a/)« e^' + . . . . 
or 

E' = A^ (aV + a\' +....)+ ^'« (a V + < V + • • • 1 f S5>. 

To illustrate by a very simple example, let 

and suppose e = 0.1 ; then, to find the mean error E of 

Y=X+X' 

we cannot take ^= i/[(0.2)2 + (O.S)^] as we should if X and X' 
were independent, but by the above formula we must take 

E = v^[(0.2)» + (0.3)« + 2 X 2 X 3 X (0.1)'] = 0.6 

as in fact we find directly, in this simple case, by first substi- 
tuting in Y the values of X and X' 
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WEIGHT OF 0B8EKVATI0KS. 

24. Observations of the same kind are said to have the same 
^ or diflferent weight according as they have the save or different 
■ v' mean (or probable) ^rrors. We assume a priori that observations 
,^ will have the samejK^iffht when they are made under precisely 
^; the same circunistancemjfiQ\\idx\\g under this designation every 
^s^ thing that can afFedk t^e observations ; but whether this condi- 
> tion has in any cage b^en realized can only be learned, a pos- 
"yj teriori^ from the meftn errors revealed by the observations them- 
selves. ■• * / \ 

In o^jder to obtain a numerical expression of the weight, let 
us suppose all our o ti e w vatioks to be compared with a standard 
fictitious oteecvafion ^e n^an error of which is any assumed 
quantity e^. Let the aqjtj^^Hfervations be subject to the mean 
error e. Let it require, a nu^ver p Xjfstamiard observations to 
-b^ combined. in-order to reoRe thA giea^ fm)r of their arith- 
metical mean to thaTofrfgiJ^u^ observation,5that is, to e ; or, 
according to (24), let 

W = e,> (36) 

then one of our actual observations is as good, that is, has the 
same weight, as p standard observations, and the number p may 
be used to denote that weight. K, in like manner, other obser- 
vations of the same kind are subject to the mean error c', and 
we have 

one of these observations has the weight of p' standard observa- 
tions, and the weights of the observations of the two actual 
series may be compared by means of the numbers p and p'. 
The weight of the fictitious observation is here the unit of 
weight; but this unit is altogether arbitrary, since it is only the 
relative weights of actual determinations that are to be con- 
sidered. 
It follows immediately, since we have 

Cj« = pe* = p'e!* 
or 

P = 7 (8^ 
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that the weights of two observatio?is are reciprocally proportional to the 
squares of their mean errors. 

The moaBuro of precision (Art. 10) and the weight are to be 
distinguished jfrom each other: the former varies inversely as 
the mean error, the latter inversely as the sq[uare of this error. 

25. To find the most probable mean of animber of obseroaliorja of 
different weights. — Let n', n", n'" .... be the given observed 
values; p', p", p'" .... their reepeetive weights. By the pre- 
ceding definition of the weight, the quantity n' may be considered 
as the mean of ^' observations of the weight unity, n" as the 
mean of p" observatioua of the weight unity, &c. We may, 
therefore, conceive the given series of observed quantities re- 
solved into a series of standard observations, all of equal weight, 
and then apply to the latter series &e principle of the arithme- 
tical mean. The whole number of equivalent standard observa- 
tions will be p' -\- p" -\- p'" -\- . ,..; the sum of the p' standard 
observations will be p'n'\ the sum of the p" standard observa^ 

K\ will be p" n" , kc. : hence the desired mean x^ will be 
ac 



'' + P 



or, more briefly, 



" + P"'»"' + -- 

hi'" + i'"'+.-- 






!*) 



This formula shows that although the above demonstration 
implies that p', p", p"' .... are whole numbers, yet any numbers, 
whole or fractional, may be used which are in the same proper. 
tion ; for / being any arbitrary factor, whole or fractional, we 
may write for (38) the following : 



^ puty V 



^ fp'n' + /p"»" + /?'"»'" + ■ 
//+//' + //"+■■.■ 



and thenfp',fp",fp"'. . . . may be regarded as the weights. 

The value of Zg is here an arithmetical mean only in the con- 
ventional sense implied in tiie substitution of fictitious observa- 
tions with uniform weights for the given observations. It may 
be called the general mean or the probable iruan. 

The weight of this general mean, referred to the unit of p',, 
j>", is = ^' + p" + p'" + 
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The mean error of the general mean will be expressed by 

where e^ is the mean error corresponding to the unit of weight. 

If Si is not given, we shall have to find it from the observations 
tliemselves. Taking the difference between Xq and each of the 
given quantities, we have the residuals 

v'=n' — Xo, v" = n" — Xq, t/" = n'" — x^y ... 

If fi', fi", €'" . . . are respectively the mean errors of w', n", n'", .... 
we shall have, as in Art. 17, 

whence 

j)'e'* = ej« = p'v^v' + y eo* 

and, in like manner, 

&c. 

The number of given values n', n" . . . being = w, the sum of 
these equations is 

?nej* = Ipw"] + IP] e,' 

which combined with the above value of e^ gives 

and consequently, also, 

ll_Spvv]_\ (4^ 

• \\(m-l)[p]/ ^ ^ 

Example. — ^Let us suppose that the observations of Saturn's 
ring in Art. 17 had been given as in the following table, where 
the mean of the first seven observations of Art. 17 is given 
= 39".179 with the weight = 7, the mean of the next following 
four = 39".285 with the weight = 4, &c. 
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p 


" 




™ 


p«v 




7 


39".179 


-0".129 


.016641 


.1165 




4 


.285 


— .023 


629 


21 




6 


.294 


-0 .014 


196 


10 




4 


.407 


+ .099 


9801 


392 




1 


.410 


+ .102 


10404 


104 




3 


.320 


+ .012 


144 


4 




3 


.377 


4-0 .069 


4761 


143 




4 


.310 


+ .002 


4 







3 


.127 


-0 .181 


32761 


983 


M = 


6 
40 


.448 
,r.= 39 .308 


+ .140 


19600 


1176 
[yi,»] = .3998 



Here the general mean a:^ found by (38) of course agrees with 
that found before. For the mean error corresponding to the 
unit of weight (which in this case is that of an observation as 
in Art. 17), we have, by (39), since m — 10, 



unit o 



and for the mean error of x„, by (40), 
_ // .3998 \_ 






■fj.. 



hi 



0".033 



which agree sufficiently weil with the foi-mer values. A perfect 
agreement in the mean errors is not to be expected, since our 
tbrmulie are based upon the supposition that we have taken a 
Bufflfieiit number of observations to exhibit the several errors 
to which they are subject in the proportion of their respective 
probabilities; and this would require a very large number of 
observations. 

26. In the application of the preceding formulse, it must be 
observed tliat when the weights of different determinations of 
the same quantity are inferred from their mean errors, we must 
be certain that there are no constant errors (that is, constant 
during the observations which compose a single determination) 
before we can combine them together according to these weights, 
unless the constant erroi-s are known to affect all the deterniina- 
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tions equally and with the same sign. For example, if ten 
measures of* the zenith distance of a star are made at one cul- 
mination, giving a mean error of 0".4, and five measures at 
another, giving a mean error of 0".8, the weights according to 
these errors would be as 4 to 1. But if it is known that the 
errors peculiar to a culmination (and aflfecting equally all the indi- 
vidual observations at that culmination) exceed 1", it would be 
better to regard the observations as of the same weight, since 
there would be a greater probability of eliminating such peculiar 
errors by taking the simple arithmetical mean. If, however, the 
observer, from considerations independent of the observations, 
can estimate the weight of determinations made under different 
circumstances, then it is evident that these weights will serve 
for the combination, if the mean accidental errors of the several 
determinations are sensibly equal. 

But if from the different circumstances we have deduced 
weights for the several determinations, and at the same time the 
mean errors (deduced from a discussion of the discrepancies of 
the observations composing each determination) are widely dif- 
ferent, it is not easy to assign any general rule for reducing the 
weights which shall not be subject to some exceptions. In such 
cases, practical observers and computers have resorted to em- 
pirical formulae, involving some arbitrary considerations, more or 
less plausible. 

In many cases we can proceed satisfactorily as follows. Let 

e = the mean accidental error of a single observation, 

7j = the mean error peculiar to a determination which rests 

upon m such observations, 
e = the total mean error of such a determination, 

then, e and tj being supposed to be independent, we shall have 

^ = ~ + ^' (41) 

m 

If then 7j can be obtained from independent considerations, this 
formula will give the value of e, and, consequently, the weight 
for each determination, and the combination may then be made 
by (38). For an example of a discussion according to these 
principles, see Vol. L Art. 236. 
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INDIRECT OBSERVATIONS, 

27. I proceed now to the application of the method of least 
squares to the solution of the general problem of determining 
the most prohable values of any number of unknown quantities 

I of which the obaen'ed quantities iire functions. The observa^ 
tioua are then said to be indirect. The particular case of direct 
observations, already considered, is, however, included in this 
general problem; being the case iu which the number of un- 
known quantities is reduced to one, and this one is directly 
observed. 

The general problem embraces two classes of problems, which 
must be distinguished from each other. In the first class, the 
unknown quantities are indepejideni, in the sense that they are 
subject to no conditions except those established by the observa- 
tions: 80 that, before taking the observaHons, any assumed system 
of values of these quantities lias the same probability as any 
other system. In the second class, there are assigned, a priori, 
certain conditions which the unknown quantities must satisfy at the 
same time that they satisfy (as nearly as possible) the conditions 
established by the ob8er\'ations. Thus, for example, if the three 
angles of a plane triangle are to be determined from observations 
of any kind, we have, a priori, the condition that the sum of 
these angles must be equal to two right angles, and all the 
systems of values which do not satisty this condition are excluded 
at the outset. This class will be briefly considered hereafter, 
under the head of *' conditioned observations ;" but our attention 
will be chiefly directed to the first class, which includes most of 
the problems occurring in astronomical inquiries. 

Again, the equations which the observations are to satisfy may 
be linear or non-linear ; the observed quantities may be explicit or 
implicit functions of the required quantities; but, for simplicity, 
we consider first the case of linear equations, to which all the 

^^^^^rs may always be reduced. 

" 28. Let us suppose the equations between the known and 

unknown quantities are of the fomi 

ax + by + c2+ + l^V 



EQUATIONS OF CONDITION FROM LINEAR FUNCTIOKS. 
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in which a, b, Cy. ...I are known quantities given by theory for 
each observation, V is the quantity observed, and x, y, 2 .... are 
the quantities to be determined. For each observation, we have 
a similar equation, and thus a system such as the following : 

a'x + 6'y +&2 '+:": + r =F' 

a"a; + fc"y + ^"2: + + V =r" 

a'"x + 6'"y + (f"2 + + r" = V"' I (^^) 

&c. &e. 

the number of these equations being greater than that of the 
unknown quantities (Art. 6). K our observations were perfect, 
all these equations would be satisfied by the same system of 

values ofXjj/yZ...; but, being imperfect, let J!f' , J!f ", J!f' " 

denote the values obtained by observation for V, F", V" 

When these values are substituted in the second members of (42), 
there will, in general, be no system of values of x, y, 2 . . . . which 
satisfies all the equations at the same time, and we can only 
determine that system which is rendered most probable by the 
observations. Let us therefore denote by iV', iV", iV'" .... the 
values which the first members of our equations obtain when 
any hypothetical or assumed system of values of x, y, 2: .... is 
substituted in them ; and put 

then 1?', v", v'" .... are the errors of the observations according 
to this hypothesis. Finally, let us put 

n'=V— M', w" = r — Jf ", n'" = T" — Jf '", .... 

then our equations may be thus expressed : 

a'x + ^'^ + (/z +^...4-/1' =1/ 
a"x + b''y +(f'2 +.'...+ n" = t/' 
a'"a: + b'"y + d"z + ....+ n'" = r'" \ (^^) 

&c. &c. 

K our observations were perfect, we should be able to find 

values of Xj T/j z . . . . which would reduce all the quantities v', t-", 

v'" .... to zero. It is usual, therefore, to write zero in the second 

members : 

a'x +b'y +</2 +....+ n' =0 

a''x + 6'V + (f'z +.... + n" = 

a'"x + 6'"y + (f''z + .... + w'" = { (^^*) 

&c. &c. 
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and these are called the equations of condition^ since they express 
the conditions which the unknown quantities are required to 
satisfy as nearly as possible. We may, however, with more rigor 
regard (43) as our equations of condition, and treat them as 
expressing the general condition that the unknown quantities 
shall be such as to give the most probable system of errors 

v', v'', I?"' 

Now, according to Art. 11, the most probable system of values 

oi x^y^ z (and, consequently, the most probable system of 

errors) is that which makes the sum of the squares of the errors 
a minimum : thus, we are to reduce to a minimum the function 

[rv] = v*v' + v"v*' + i;"V" + . . . . 

Regarding [?w] as a function of the variables x^y^z... (which we 
must remember are here independent), the condition of minimum 
requires that its derivatives taken with reference to each variable 
shall each be zero ; that is, 

^CHHI^q IM^Q ^M^Q _ 

dx ^ dy ^ dz 

or 

dx dx dx 

dy ^ dy ^ dy ^ } (44) 

az az az 

&c. 

(which we might have obtained directly from (10) by substituting 
y)' J == kJ = kv, and dividing by the constant k). But, by diflfer- 
entiating the equations (43) with reference to Xy y^ z . . . . succes- 
sively, we have 

dv' , dv' ,, dv' . 

dx ' dy ' dz ' '' 

dv" _ „ dv" _ „ ^^"_v^ 

-- — a y — - — , -- — cr , . ..• 

dx dy dz 

&c. &c. &c. 

BO that (44) are the same as the following : 
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ft't,' _|. 6'V' + fc'"i7'" + .... = 
drf + d*v" + c"V" + .... = ( C*^*) 

&c. 

The number of these equations is the same as that of the un- 
known quantities ; and if we now substitute in them the, values 
of I?', v", v'" . . . from (43), we have the final or, as we shall call 
them, the normal equations, which determine the most probable 
values of X, 3/, 2 ... . 

NORMAL EQUATIONS. 

29. We see by (44*) that to form the first normal equation we 
multiply each of the equations of condition (43) or (43*) by the 
coeflicient of x in that equation, and then form the sum of all 
the equations thus multiplied. The resulting equation is called 
the normal equation in x.* The sum of the equations of condi- 
tion severally multiplied by the coefficients of y is the normal 
equation in y, &c. To abbreviate the expression of these sums, 
we put 

[aa] = a' a' + a'' a'' + *a'"a'" + . . . . 

lab-] = a'y + a"6" + a'"6'" + . . . . 

\ac\ = aid + a'V + a"V" + . . . . 

&c. &c. 

then the normal equations are 



\aa\ X + [a6] y + \ac\ z -^ — -j- [an] = 
[a6] X + [66] y + [6c] ^: + . . . . + [6n] = 

\ac\ X + [6c] y + \cc\ z -\- -f [<^] = ^ 

&c. &c. 



(45) 



80. The formation of such normal equations is one of tiie 
most laborious parts of the computations involved in the method 
of least squares, especially when the number of equations is very 
great. It is important to have a means of verification, or 
"control,*' to insure their accuracy, before proceeding with the 
next important process of elimination. A very simple and 
effective control is the following. 



^ The '* normal equation in a;" is so called because it is the equation which deter- 
mines the most probable value of x when the other variables are reduced to zero, <ff 
when X is the only unknown quantity ; and so of the others. 
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Form the sums of the coefficients of the unknown quantities 
in the several equations, namely, 



b 



a' +6' + c' -H . 
a" + b" + <;'■ + . 
«'" + i'" + C" + ■ 



C4fi) 



If we multiply each of these by its n, and add the prntlucta, we 
have 

Also, multiplying each of (46) by its a, and adding, then each 
by it s b, and adding, and ao on, we have 

[aa] + lab-] + [cc] + . . . . = [<7s] \ 

[ah-] +[66] +[6c] +.... = [*.] 

[„c] + [6c] +[cc] +.,.. = [cO C-IS) 



The equations (47) must be satisfied when the absolute terms of 
the normal equations are correct, and (48) when the coefficients 
of the unknown quantities are correct. 

31. The normal equations will give determinate values of 

3; y, z , provided they are really independent. If, however. 

any two of them become identical by the multiplication of either 
of them by a constant, the number of independent equations is, 
ill fact, one less than that of the unknown quantities, and the 
problem becomes indeterminate. This difficulty does not arise 
from the method by which the normal equations are formed, but 
from the nature of the given equations of condition. In any 
such case, additional observations are necessary, for which the 
coefficients have such varied values as to lead to independent 
equations. Even when t^vo equations cannot be reduced pre- 
cisely to a single one by the introduction of a constant factor, if 
tliey can be made very nearly identical, the problem is still pi-ac- 
tically indeterminate. The indetermi nation will become evident 
in the actual elimination in practice when any one of the un- 
known quantities comes out with ao small a coefficient that small 
errors in the observations would greatly change this coefficient, 
CSee Art. 52.) 
Vol. II.— 33 
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32. By whatever method the elimination is performed, we 
shall necessarily arrive at the same final values of the unknown 
quantities ; but, when the number of equations is considerable, 
the method of substitution, with Gauss's convenient notation, is 
universally followed; but, for the present, leaving the reader to 
choose his method, I proceed to explain the principles by which 
the mean errors of the values of x, y^ z . , , . are determined. 

MEAN ERRORS AND WEIGHTS OF THE UNKNOWN QUANTITIES. 

33. Since we have put n' = V —M\ ri" = U'—M'', &c. (Art. 28), 
the mean error of n\ n", w'". ... is also that of M% J!f ", Jtf' ", . . . . ; 

that is, the mean error of n', n", n'" is to be regarded as the 

mean error of an observation. If the elimination of the normal 
equations were fully carried out, each unknown quantity would 

be finally expressed as a linear function of 7i', r?/', w"', , and the 

mean errors of the latter being given, those of the unknown 
quantities would follow by the principle of Art. 20. It results, 
however, from the symmetry of the normal equations that several 
forms may be obtained for computing directly the weights of the 
unknown quantities, and from these weights the mean errors 
can afterwards be found. 

34. First method of computing the weights of the unknown quantities. 

— For simplicity, let us first suppose all the observations to be 

of equal weight, or the mean errors of n', n", n'" to be equal. 

Let 

e = the mean error of an observation, 

e = the mean «rror of the value of x found from the normal 

SB 

equations, 
jp^= the weight of the value of Xy the weight of an observa- 
tion being unity; 

then (Art. 24) 



X 



Now, let us suppose the elimination to be performed by the 
method of indeterminate coeflicients. Let the first equation of 
(45) be multiplied by Q^ the second by §', the third by §", &c., 
and the products added. Then let the factors §, ^, Q" .... 
(whose number is the same as that of the unknown quantities) 
be supposed to be determined so that in this final equation the 
coefiicients of all the unknown quantities shall be zero, except 



■^1 



(aO^. ; U ^ C^ 'r^ J K\ - ■■::' ^^' -^ 
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that of x^ which shall be unity. The conditions for determining 
these factors are, therefore, 

\aa\ C + \ph'\ ^ -f [rtc] ^' + ....= 1 
[a6] e + [66 ] Q' + [6c] ^' + .... = 
[ac] C + [6c] 6' + [cc] §" + . . . . = 
&c. &c. 



(49) 






and the final equation in x is 

X + [a?i] § + [6n] Q' + [en] §" + .... = (50) 

Comparing (45) and (49), we see that the coefficients of 

§, Q', §".... are the same as those of x^y^z , but that the 

absolute terms are — 1 in (49) instead of [an] in (45), and zero 
instead of [6n], [en], &c. Hence, if the elimination of (45) were 

carried out, and the values oi x^y^ z determined in terms of 

n', n", n'". ..., the values of §, Q', §" would be found from 

these by merely putting [an\ = — 1, and [6n] = [en], &c. = 0. 
This is also evident from (50). I shall now show that Q is the 
reciprocal of the required weight of x. 

The final value of x being a linear function of n', n", n'". . . ., 
the equation (50) may be supposed to be developed in the form 

X -f aW + a"n" + a'"w'" + ....= (51) 

in which a', a", a'", are functions of a', 6', , a", 6", , &c. ; 

and these functions are immediately found by developing [an], 
[An], &c., in (50) ; for we then have, by comparing the coefficients 
of (50) and (51), 

a' = a' Q + b' Q' + e' §" + .... 

a'"= a'" Q + 6'" §' + e"' §" + ... . ^ ^ ^ 

&c. &c. 

Multiplying each of these equations by its a, and adding all the 
products, we obtain, by (49), 

a a -\- d a -f- O, a -\- . . . , = 1 

Multiplying each of (52) by its 6, and adding, we obtain, by (49), 

6'tt' + 6"a" + 6'"a'" +....= 

and so on for as many equations as there are unknown quantities. 
These relations are briefly expressed thus: 

[aa] = 1 [6a] = [ca] = 0, &c. (53) 
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If, then, each of (52) is multiplied by its a, and the results are 
added, we find, by (53), 

[aa] = a''+ a"*+ *'"'+....= Q (64) 

But, by Art. 20, when e is the mean error of each of the quan- 
tities n', n", n'", . . . ., the mean error of x found by (51) is 

e,= e/[aa] 

Ilence 



e» 1 1 



as was to be proved. 

Hence we have a first method of finding the weights. In the 
first normal equation write — 1 for the absolute term [flf??], arul in the 
other equations zero for each of the absolute terms [6w], [«i], &c. ; the 
value of X then found from these equations will be the reciprocal of the 
weight of the value of x found by the general elimination. 

This rule is to be applied to each of the unknown quantities 
in succession, so that the reciprocal of the weight of y is that 
value lOf.y which will be found by putting '^n\-= — 1, and 
[an] = [en] = &c. = ; the reciprocal of the weight of z is that 
value of z which will be found by putting [en] = — 1, and 
[an] = [6n], &c. = ; &c. 

It is evident, moreover, that although we have deduced the 
rule by the use of indeterminate multipliers, it must hold good 
whatever method of elimination is adopted. 

35. Second method of computing the weights of the unknown quan- 
tities, — ^K we write the normal equations thus, 

[aa] X + [P'b'] y + [aC] z + — + [^w] = ^ 
[a6] x + [66] y + U>c]z +....+ [6w] = B 
[ac] X + [6e] y + [cc] z + — + [en] = C 
&c. &c. 

and perform the elimination, we shall obtain a:,y, 2'.... in terms 
of [aa], [a6], &c., and of A, jB, (7, &c. ; and if in the general values 
thus found we make ^ = jB = C, &c. = 0, these values will be 
reduced to those which would be found by carrying out the 
elimination with zero in the second .members of the normal 
equations. If we suppose the elimination performed by means 



\ 
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of the indeterminate factors Q, Q', Q" already employed, the 

final equation for determining x will be 

a: + [«"]« + [H9'+ Me"+ ■■■• =6'! +§'-8 +e"c+.... 

where the coefficient of A is the reciprocal of the required weight 
of X, But, whatever method of elimination is employed, the 
coefficient of A in this general value of x will neceeearily be the 
same; and hence we derive the eeeond method of determining 
the weights : Write A, B, C, &c., iTistead of 0, in the second members 
of Ike normal equations, and rarry out the elinmaiion (by any method 

at pleasure) ; then the filial values oj x,y, z are those terms in the 

general values which are independent of A, B, C ; the weight of x 

is the reeiprocal of the coefficient of A in the general value of x; the 
weight of y is the reciprocal of the coefficient of B in the general value 
ofy; &c. 

36. Third method of computing the weights of the unknown quantities. 
— Let U8 suppose the elimination to be performed by the method 

of substitution, still reta.ining A,B, C in the second members, 

as in the preceding article. The final equation in x, accoi'ding 
to this method, is found by substituting in the first normal equa- 
tion the values ofy, 2.,.. given by the other equations. These 
substitutions do not affect the coefficient of A, which remains 
unity, so long as no reduction is made after the substitutions. 
Thus, the final equation in x is of the form 

Rr= T+ A + tcrma in 5, C, ... . 

which T is the sum of all the absolute quantities resulting 

from the substitution, and is a function of [«(*], [afi], [«ji]. 

Hence the value of x ia 



^Thu 

froE 
HeL 



in which ^ is the final value of x which results when A — B 
^ C . ■ . = 0, and -^ is necessarily the quantity denoted by Q in 

the preceding articles. Therefore E is the weight of x, and 
hence we have a third method of finding the weights : Let the first 
mJiorinal equation (the equation in x, Art. 29) be taken as the fined 
equation for determining x, and substitute in it the values of y, z in 
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terms of x as found from the, remaining equations; then, before freeing 
the equation of fractions or introducing any reduction factor , the coeffi- 
cient of X in this equation is the weight of the value of x. In the same 
manner, substitute in the second normal equation (the equation in y) 
the values of x,z,...in terms of y as found from the other equa- 
tions; the coefficient of y is then the weight of the value of y; and so 
proceed for each unknown quantity. 

According to this method we determine each unknown quan- 
tity, together with its weight, by a separate elimination carried 
through all the equations, in each case changing the order of 
elimination, until every unknown quantity has been made to 
come out the last. The algorithm of this process, with Gauss's 
convenient system of notation, will be given hereafter (Art. 45). 

37. To find the mean error of observation. — The weight of x being 
found, we have the ratio of e^ to e, but we have yet to determine 
e, which, in general, cannot be assigned a priori, but must be 
deduced a posteriori, that is, from the observations, and conse- 
quently from the equations of condition. The residuals v',i?",t?'"...., 
in (43), are those which result when the most probable values of 
x,y,z.... (namely, those resulting from the normal equations) 
are substituted in the first members. The actual or true errors 
(Art. 17) of observation are, however, those valuea of the first 
members of (43) which result when the true values ofx,y,z,.... 
are substituted. 

Let X + AX, y + ^y, z + ^z,. . . .he the true values which, sub- 
stituted in the equations of condition, give the true residuals 
u% m", m'". . . . ; so that we have 

a' (x 4- ax) + b^ (y + Ay) + (/ (z + az) + n' =u' 

a" Ix + AX) + b" ly + Ay) + c" (2r + a-j) + . . . -. %" = u" . 
a'"(^ + ax) + U\y + Ay) + d"{z + a^) + . . . . n'"= v!"\ (^^) 

&c. &c. 

If these equations be multiplied by a', a", a'". . . ., respectively, 
the sum of the products is 

\aa'\x +[«*3y +M^ +....+ [an] l^^r^^-i 
-}- \aa'\ AX -j- [a6] Ay + [ac] az + J ^ ■* 

which by the first of (45) is reduced to 

[aa"] AX -f [a6] Ay -f [ac] a^ + — — [au] = 
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In the same manner, multiplying each of the equations (56) by 
its by c, &c., successively, we form the other equations of the 
following group : 

[ad] A.r + [^^] ^y + [^^] ^^ + — [^w] = 

[a6] Aic + [bb^ /ay + [ftc] a;? + . . . . — [ftw] = . 

[ac] AX -f [ftc] Ay + [cc] A2 + — [ct^] = ' ^ "^ 

&e. &e. 

These being of the same form as the normal equations (45), we 
see that the value of ax resulting from them will be of the samq 
form as that of x resulting from (45), with only the substitution 
of — u for n: hence, by (51), 

AX — a'w' — a"M" — a"V" — ....= (58) 

Again, multiplying (56) by v\ v", v'". , . ., respectively, the sum 
of the products is, by (44*), reduced to 

[vn] = [vm] 

and in the same manner, from (43), 

[vn] = [vu] 
whence 

[vm] = [vv] = [vri] (59) 

The sum of the products obtained by multiplying the equations 
(43) respectively by u\ w", w"'. ... is 

[au] X + [6m] y + [cu] z -{■ + \nu] = [vu] = [w] 

and from (56), in the same manner, 

{au]x +[bu\y + [cw] 2r +..•• + C^w] ) ., 
+ [«w] AX + \bu\ Ay + [cm] a^ + j ^ 

which two equations give 

[mm] = [tn?] + [au] AX -|- \bu] Ay + [cm] A2r + (60) 

Now, \mi\ being the sum of the squares of the true errors of the 
observations, its value is, as in Art. 17, = mse, if we put 

m = the number of observations, 
•= the number of equations of condition. 
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Consequently, if we could assume ax, a^ .... to vanish, we should 

have 

[tw] 

€€ = = 

m 

and this will usually give a close approximation to the value of 
e, but it will give the true value only in the exceedingly impro- 
bable case in which the values of a:, y, 2:. . . . are absolutely true, 
whereas they are to be regarded only as the most probable ones 
furnished by the observations. This formula, then, must always 
give too small a value of e, since it ascribes too high a degree of 
precision to the observations. We must, therefore, add to [rr] 
the quantities [au] ao:, [6u] a^, &c., as in (60) ; but, as we cannot 
assign any other than approximate values of these quantities, let 
us assume for them their mean values as found by the theory of 
mean errors. The mean value of [au] ax will be found by mul- 
tiplying together 

[au] = a'u' + a"u" + a'"M'" + . . . . 

and AX = a'w' + a"M" + a'"w'" + . . . . 

observing that the errors m', m", w'". . . ., when we consider only 
their mean values, are to be regarded as having the double sign 
±: ; 80 that the mean value of the product will contain only the 
terms a' a' m' u', a" a'' w" m", &c. Hence we take 



(XUj AX = (laUU'-\-(laUU -\- CI a U U -f-.... 

and substituting in this the mean value of u^u\ u''u'\ fcc, which 
in each case is ££, we have 

[au\ AX = (a! a! + a" a." + a!" a!" + ....)«« 

or, finally, by (63), 

[au] AX = ee 

In the same manner, it must follow that ee is the mean value of 
each of the terms [bu] Ay, [cii] A^, &c. If then we put 

/I = the number of unknown quantities, 

the equation (60) becomes 



METHOD OP LEAST SQUABE8. 



[""] 



(61) 



It is to be observed that when there is but one unknown 
quantity, or /^ = 1, this general form is reduced to the simple 
one (26), already given for direct observations. 

Finally, p^, p^, p^, denoting the weights of x, ^, 2.... found 

by any of the preceding methods, we have 



Vp, 



, 4o. 



38. 



Example. — Let ua snppoee the following veiy e; 
of condition to be given ;* 



. (62) 
nple equa- 



3x + 2y — 5z — 5 ^ 

43; -|- 1/ + 4i — 21 = 

— x 4- 8y -f- 33 — 14 ^ 







If but the first three of these equations had been given, the 
problem would have been determinate, "We should find from 

18 23 13 

them a; = — -,y — — , z — — -, and we should have to accept these 

values as final ones, with no means of judging of their accuracy, 
or of that of the observations upon which the equations are sup- 
posed to depend. A fourth observation having given ua our 
fourth equation, we find that the values of x,y, z derived from 
the first three will not satisfy it, for when they are substituted in 
it the first member becomes — ^, instead of zero. If we deter- 
mine the values of z, ^, and z from any three of the equations, 
and substitute these values in the fourth, we shall find a residuah 
Each one of the four systems of values of the unknown quantities 
thus found satisfies three equations exactly, and the fourth 
approximately; but, all the observations being subject to error, 
the most probable system of values can seldom satisfy any one 
of the equations exactly. Hence the necessity of a principle of 
computation which shall lead as directly as possible to such a 
probable system of values; and this principle is furnished by the 
pjethod of least squares. 



• Gausb, Theo 
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We are, then, by Art. 29, to deduce from these four equations 
three normal equations, and the values of a:, y, z which exactly 
satisfy these are to be regarded as the most probable values. 

To form the first normal equation, we multiply the first of the 
above equations of condition by 1 (= a'), the second by 3 (= a"), 
the third by 4 (=a'"), and the fourth by — 1 (=«*0, and add the 
products. We thus find [oa] = 27, [a6] = 6, \ac\ = 0, and 
\^m-\ = — 88. 

To form the second normal equation, we multiply the first 
equation of condition by — 1 (= ^0) ^^® second by 2 (= 6"), the 
third by 1 (= V\ and the fourth by 3 (= U% and add the pro- 
ducts. Ve thus find [a6] = 6, [66] = 15, \hc\ = 1, \hn\ = — 70. 

The third normal equation is formed by multiplying the first 
equation of condition by 2 (= c'), the second by — 5 (= c")> *^^ 
third by 4 (== c'''), and the fourth by 3 (= c*^), and adding the 
products. We find [ac] = 0, [be] = 1, [cc'] = 54, [cri] = — 107. 

Hence our normal equations are 

27x + 6y — 88 = 

6x + lby + z— 70 = 

y _)_ 542: _ 107 = 

the solution of which gives, as the most probable values, 

49154 ^._ 

X = = 2.470 

19899 

y = ^?lL^3.551 
^ 737 

6633 

In order to determine the mean, and hence also the probable, 
errors of these values, let us first determine their weights accord- 
ing to the preceding methods. 

First. By the method of Art. 34, we first write — 1, 0, , for 
the absolute terms of the three normal equations, and we have 
the three equations for determining the weight of x, 

27a/ + 0}/— 1 = 

6x!-\-lb\/+ zf = 

y'+54/ = 

in which accents are employed to distinguish the particular 
values from the above general ones. These give 
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^^ 809 

19899 
which is the reciprocal of the required weight. Hence, 

» = = 24.597 

^' 809 

In a similar manner, to find the weight of y, we take the 
equations 

27a:" + 6/' =0 

6a/'+15y"+ 2"— 1 = 

y" + 542" =0 



and find 



whence 



y = — 

^ 737 



P.= fj = ^^-^^ 



And to find the weight of ^, the equations 

27a/" + Gy"' =0 

Gx"' + 15y'" + y" =0 
y'"+542'"— 1 = 
which give 

y"=ii- 
2211 

and 

«=??11 = 53.927 
^' 41 

Secondly. By the method of Art. 35, we write our normal 

equations thus : 

27a: + 6y — 88=^ 

6a: 4- 15y + z— ^^ = B 
y 4- 542: — 107 = a 

and, carrying out the elimination as if J., Ji5, and C7 were known 
quantities, we find 

lj»899 X = 49154 + (809)A — 324 5 + 6 

. '• 737y = 2617 — 12 A + (54)5 — G 

6633^ = 12707 + 2il — 9 -B + (123) C 



524 APPEITDIX. 

and, therefore, 

49154 . . ,, . , , 19899 

X = with the weight o = -— -- 

19899 * '^" 809 

2617 ,, ,, ,, 787 

^ 737 ^^ 64 

12707 ,, ,, ,, 6633 

6633 -^' 123 

the same as by the first method. 

Thirdly. By the method of Art. 36, to find x and its weight 
we eliminate y and z from the equation in x (the first normal 
equation) by means of the other equations, employing successive 
substitutions. The last normal equation gives 

1 , 107 

z = y A 

54 ^ ^ 54 

which being substituted in the second gives 

^ ,809 3673 . 
%x A y = 

^ 54 ^ 54 
The value of y from this, namely, 

324 . 3673 

v = x 4- 

^ 809 ^ 809 

being substituted in the first normal equation, and no reduction 

being made, gives 

19899 49154 _ 

809 ^ 809 "■ 

where the coefficient of x is the weight, and the value of x is the 
same as before found. 

To find y and its weight, we make the second the final equa- 
tion. From the first and third we find 

— _-^ 88 

^ — 27 ^ "•" 27 



^ = -£iy + 



107 



54^ ' 54 
which substituted in the second give 

737 2617 ^ 

■54"y--5r = ^ 

where the coefficient of y is its weight. 
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Finally, to find z lyith its weight, we make the third ziormal 
equation the final one. From the first two we find 



y = — 



_9_ 454 

123 ^ "^ 128 



which substituted in the third gives 



6633 ^ 12707 
z — 



128 



123 







where the coeflicient of z is its weight, and its value is the same 
as was before found. 

By a little attention, it will be perceived that the three methods 
involve essentially the same numerical operations. 

We are next to find the mean errors of a:, y, and z\ for which 
purpose we must first find the giean^error of an ohaervation 






assuming here, for the sake of illustration, that the absolute terms 
of the given equations of condition are the observed quantities, 
and that they are subject to the same mean error. Substituting 
in these equations the above found values of x, y^ and z^ we 
obtain the residuals as follows : 



No. 





w 


1 


— 0.249 


0.0620 


2 


— 0.068 


.0046 


8 


+ 0.095 


.0090 


4 


— 0.069 


.0048 



Hence, by (61), 



m =-- 4, A = 8, [Dv] = 0.0804 

m — fi ' 



0.0804 



c = 1/0.0804 = 0.284 



which is the mean error of an observation, so far as this error 
can be inferred from so small a number^of observations. . ^^See 
the next article.) Consequently, the mean errors of x, tf, and z 
are as follows : 



e 



€. = -f- = 0.057 
r = -4- = 0.077 



e. 



Vp, 



0.089 
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Multiplying these errors by the constant 0.6745, we shall have 
(Art. 15) the probable errors as follows: 

Probable error of an observation >= 0.192 
« " a: = 0.038 

« « y =0.052 

" *< z =0.026 

39. It has already been remarked in the foregoing pages, and 
the remark is especially important in the present connection, 
that the method of least squares supposes in general a great 
number of observations to have been taken, or a number suffi- 
ciently great to determine approximately the errora to which the 
observations are liable. Theoretically, the greater the number 
of observations the more nearly will the series of residuals ex- 
press the series of actual errors, and, consequently, the more 
correct will be the value of e inferred from these residuals. In 
practice, therefore, no dependence should be placed upon the 
mean or probable errors deduced from so small a number of 
observations as we have employed, for the sake of brevity and 
clearness, in the preceding example. Nevertheless, the method 
is, even in this case, the best adapted for determining the most 
probable values of the unknown quantities deducible from the 
given observations, and also their relative degree of precision. 
Thus, in this example, the degrees of precision (denoted by A, 
Art. 10) of x, y, and ^, being inversely proportional to the mean 
errors, or directly proportional to the square roots of the weights, 
are nearly as the numbers 5, 3.7, and 7.3, so that from the four 
given observations z is about twice as accurately found as y, 
while the precision of x falls between that of y and z. But we 
can place but little dependence upon the result which assigns 
0.284 as the mean error of observation, and 0.057, 0.077, 0.039 
as the mean errors of x^ y, and e, because this result is derived 
from too small a number of observations. 

EQUATIONS OF CONDITION FROM NON-LINEAR FUNCTIONS. 

40. Let the relation between the observed quantities V, V'/, 
V". . . . and the unknown quantities X, I^ ^. . . . be, for the ob- 
servations severally, 

/' (7', X,Y,Z, ) = 

/"("T", X,Y,Z, ) = 

/'"(F'", XjY.Z, ) = ( (^) 

&c. 
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Let the values of V, V", V" . ..., fouod by obaervation, be 
M', M", M'" .... These values being substituted, we ehall 
have the equations 



/' (M', XY.Z, ■■■■') =0 

f" (if", j:,y,z,....) = q 

f"'{M"',X,Y,Z,....-)=Q 



(64) 



from which the values of X, Y, Z are to be found. But, aa 

we c-auiiot effect the direct solution of these equations according 
to the method of least squares ao long as they are not linear, we 
resort to the following indirect process, by which linear equations 

of condition are formed. Let fl/ysroxiwiofe values of X, Y, Z be 

found, either by some independent method or from a sufficient 
number of the equations (64) treated by any suitable process, and 

denote these approximate values by Xq, Y^, Z^ Let the most 

probable values be 



,r^^^ + a:, 



= Y, + y, 



Z=Z^-i-s,. 



then X, y, z,. . . are the corrections required to reduce our ap- 
proximate values to the most probable values; in other words, 

x,y,z are the most probable eorreetions of the approximate 

values, and the method of least squares is now to be applied in 
finding these corrections. 

Substitute the approximate values Xo, Y^, Z^. . .in (63), and 
find, by resoh-ing the equations, the corresponding values of 
V, V" .... which denote by F/, VJ," .... These will be func- 
tions which may be thus generally expressed: 

f; = F' (X„, y„, 2o . . . .) 
&c. 

Now, the values of T', V" .... which result when the most pro- 
bable values Xb + 37, Yf,-\- y, Za-^ z are substituted, and which 
are yet unknown, being denoted by N', N" .... we have 



N' = F' (A', + .r, 
if"^ir"(j; + x, 



Y„+y,Z„ + z,....) 



and by Taylor's Theorem, when wc neglect the higher powers 
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of a:, y, 2r . . . . which arc supposed to be very small quantities, we 
have 

&c. &c. 

dV dV" dV* dV" 

where -T-sf' ——-» &c., -r^ ttt"' &c. are simply the values of the 
dX^ dX^ ' dY^ dY^ ^ ^ 

derivatives of F', V" .... found by differentiating (68) with 
reference to each of the variables, and afterwards substituting 
X^^ !Fi, &c. for X, F, . . . . &c. 

K now we denote the derivatives of Y\ F" .... with reference 
to Xby a', a" .... ; their derivatives with reference to Fby 6', 
6'' . . . . &c. : so that 

W =Yl + rt'a; H- 6'y + (/^ + . . . . 

&e. &c. 

and then also put 

'^ = W-^ M', v" =N'' — M", &c. 

n' = Fo' — M', n" = Fo" — Jf ", &c. 

our equations become 

a!x + 6'y + dz + . . . + n' :±± v' ■ • 
fr"a? + ft"y + c"2: + . . . + n" == t?" 
a'"x + 6'"y + 0"'^ + . . . + n'" = i/" 
&c. &c. 

in which «', 6' . . . a", 6" . . . n', n" . . . are all known quantities ; 
and v\ v" .... are the residual errors of observation. These 
equations of condition are precisely like those already treated, 
and, being solved by the same method, give the most probable 
values of a:, y, ^ . . . ., and hence, also, the most probable values 
ofJ',F,Z.... 

' This process rests upon the assumptiom that the approximate 
values Jf^, Fo, Z^... are already so nearly correct that the squares 
oi Xy y^ z . . . may be neglected. But should the values found 
for x^y^z . . . show that this assumption was not admissible^ the 
' computation is to be repeated, starting with the last found values 
}[q'\- x^ yj, + y, ^ + ^ . . . as the approximate values ; and tten 
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the corrcetione which these last require will generally be ho small 

, that their liigher powers may be neglected without sensible erroi'. 

I However, should this still not be tJie caKe, successive approxima- 
tions, commencing always with the last found values, will at 

J Leiigth lead to values whicli require only corrections suitably 

^ small. 

I Even when the given function is already linear, it ia mostly 

II expedient to follow the general nnethod just given : namely, to 
I substitute approximate values and form ecjuations of condition 

to determine their corrections. This reduces t, y, z . . . to small 
quantities, greatly simplifies the computations, and diminiehea 
the chauce of error. 

TREATMENT 01 EQUATIONS OF CONDITION WHEN THB OBSERVATIONS 
HAVE DIFFERENT WEIGUTS. 

41. The process above explained assumes that all the obser\'a- 
1 tions are subject to the same mean error, and hence are all of 
I the same weight. The more general case, in which the obser- 
vations are of different weights, is easily reduced to this simple 
case. For, let 

a'x + b'y + (/£ + . . . . + n' = t^ 

be an equation of eonditiou of the weight p' ; that is, one formed 
for an observation of the weight p'. The mean error of an ob- 
servation of the weight unity being e,, the mean error of the 
actual observation, and, therefore, also ofn', is e' — — '-r Hence 
the mean error of n'y^p' is, by Art. 20, equal to e'i />', that ia, 
equal to e,. If, therefore, we multiply the equation by yp', so 
that we have 

a'l/p' . X + fe'vV- y + ^V¥- ^ + ■ ■ ■ + "V/ — "V P 
it becomes an equation in which the mean error of the absolute 
term is the mean error of an observation of the weight unity. 
( Hence we have only to multiply each equatiou of condition by 
the square root of its weight in order to reduce them all to the 
same unit of weight ; after which the normal equations will be 
found ae in other cases. 

The mean error of observation, found by (61) from the equa- 
il tions of condition thus transformed, will be that of an observa- 

Vot,. II.— 34 
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tion of the weight unity, and the weights of the unknown quan- 
tities will come out with reference to the same unit. 



ELIMINATION OF THE UNKNOWN QUANTITIES FROM THE NORMAL 
EQUATIONS BY THE METHOD OF SUBSTITUTION, ACCORDING TO 
QAUSS. 

42. By means of a peculiar notation proposed by Gauss, the 
elimination by substitution is carried on so as to preserve 
throughout the symmetry which exists in the normal equations. 
In order to explain this method, it will be expedient to suppose 
a limited number of unknown quantities. I shall take but/owr, 
but shall give the process in so general a form that it may readily 
be extended to any number. 

The unknown quantities will be denoted by 

X, y, z, Wy 

and their coefficients in the equations of condition by 

a, h^ c, dy 

respectively, with sub-numerals denoting the number of the 
equation or observation upon which it depends, and by 

^i> Wa» ^8> *^- 

the absolute terms of the 1st, 2d, 3d, &c. equations respectively : 
so that the m equations of condition (here supposed to be 
reduced to the same weight by Art. 41) will be 

«i^ + KV + o^z + rfjW? + Wj = 
a^ + ^sy + V + ^^ + w^ = 

«8^ + *ay + V + ^z^ + Ws = \ (65) 

• •• ••• ••• 

• •• ••• •«• 

a^ + ^«y + <^J^ + ^^ + w« = ^ 
and the four normal equations formed fropi these are 



[aa] X -f [a6] y -f \ac\ z + [a^ w + \an\ =0 
[a6] X + [66] y + [6c] i + \bd:\ w + [6n] = 
lac'\x + [be'] y + Ice] z + led] w -t [cwj = 



lad] X + lbd]y + led] z + Idd] w ^ Idn] = 



^Idn] = 



(66) 
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The value of x from the first equation is 

[a5] [flc] \a^^ \an\ 

[aa\ {aa\ [aa] [aa\ 

If this is substituted in the other three equations, we shall pre- 
serve the symmetry of the result by the following notation : 



\aa] 

[aa\ 

\aa] 

The three equations thus become 

\bh.V]y +\hc.r\z+ [6^. !]«/? + [6w.l] = 
[6<? . 1] y + [cc . 1] 2r -j- [cd .V]w -\- 
[6(i.l]y + [cd.l]2r + [rf(f.l]M7 + 



\aa\ 

[bri] — E^ [an] = [6w . 1] 
[aa] 

['^^ - £^ ^""^ = [en . 1] 
[aa] 



[en . 1] = I 
[dnA] = ) 



(67) 



The presence of the numeral 1 is all that distinguishes these 
from original normal equations in y, 2, and iv. The elimination 
of y will, therefore, be effected in the same manner as that of x. 
Thus, from the first, we have 

_ [be . 1] ^ [bd . 1] [bn . 1] 



y = 



[66 . 1] [66 . 1] 



[66 . 1] 



the substitution of which in the other two equations leads to the 
following notation : 



[cc.l] — 



[cdA]^ 
[dd.l]^ 



[be A] 
[bb . 1] 

[6c. 1] 
[66 . 1] 
[bd . 1] 
[66.1] 



[6c.l] = [cc.2] 
[bdA] = [cd,2] 



[bdA] = [dd.2] 



{en . 1] — 
[rfn.l] — 
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and the resalting equations are 



[cc . 2] z + [cd . 2] M? + [(?n . 2] = 
lcd,2]z +[dd,2]w + [dn.2] =0 



} (68) 



From the first of these we have 



r^ ^^-2]^ [en. 2] 
Ice . 2] Ice . 2] 

which, sabstitated in the second, leads to the following notation : 



[<W . 2] ~ i^^^ [ci . 2] = [dfrf .3] 
and the resalting equation is 



[(fn.2]— ^f^L?3 [en. 2] = [e^n.3] 



whence 



.^ 



[<W.3Ju7 + [rfn.3] = 



[/Zn . 3] 
Idd . 3] 



(69) 




Having thus foun^ |ftf^j^9|pp^stitute its value in the first of (68), 
and deduce. 2. T!ieif.tifel^ues of z and w being substituted in 
the first of {Sf^ wie deduce y ; and finally, substituting the values 
j/, Zy and w in the first of (66), we deduce z. These latter substi- 
tutions are made in the numerical computation, but it is not 
l^eceesary to write out here the formula which result from the 
litei^l substitutions, as it would not facilitate the computation. 

It may be observed that all the auxiliaries [66 . 1], [6c . 1], [cc . 2], 
&c., may be expressed by the general formula 

a, ^, Y denoting any three letters, and /z any numeral. 

For the convenience of reference, the final equations employed 
in the actual computation are brought together as follows, thje 
coefficient of that unknown quantity which is found from each 
after the substitution of the values of the others being reduced 
to unity : 
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M' 






[no] 



[fa.l] [M.l] 

"[».l] '■'"[4S.1] 



I" i«o] 



(70) 



^ [CC.2J ^[ra.2] 

Ab the number of unknown quantities increases, tlie number of 
auxiliaries to be fouud increases very rapidly. If we include the 
coefficients and absolute terms of the normal equations, the 
whole number of auxiliaiies is shown in the following scheme :* 



No. of unknown qnantities 


1 


2 


3 


4 1 5 j 6 


7 


8 




'' 


7 


.e 


30 1 50 1 77 


112 









48. For the purpose of verification, it is expedient to repeat 

the elimination in inverse order, commencing with the last 
normal equation and ending with the first, which will bring out x. 
It will not be neccseary to write out the formulte for this inverse 
elimination, since when the form for computation has been once 
prepared, it sufBces to place in it the coefficients of the normal 
equations in inverse order, and then to proceed with the numeri- 
cal operations precisely as in the first elimination. The unknown 
quantities coming out in the fijst elimination in the order w, 2, 
^, X, they will in the second come out in the order x, y, z, u\ 

This inversion has rIbo the advantage of giving the weights of 
all the unknown quantities with the greatest facility, as will 
hereafter be shown. 



44. A very complete final verification, or " control," is obtained 
as follows. Substitute the values of x, y, z, w in the equations of 
condition, and thus find the residuals v,, Uj, i\. . . . v^, or the 
values which the first members assume. Ponu the sum 



M = 



I * Tlie Dumber of a. 



lilittriea will be, in general, 

.•(■■+ !)(.■+ 6) 



■rfaere i deaotes Lhe a 
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which is also required in finding the mean error of observation 
by (61). Also form the following new auxiliaries : 



Inn} = n^w, + nj^^ + nJ^^+ + nji^ 

rnn.2]-£??i^ = [nn.3] 



[aaj 



- -_ [bn.lV - .„ 



[nn . 3] — -^ = [nn . 4] 



then, if the whole computation, both of the normal equations 
themselves and of the subsequent elimination, is correct, we 
must have 

[rr] = [nn . 4] (71) 

To demonstrate this, we observe first that we have already, by (69), 

If now we go back to the equations of condition, and multiply 
each by its n, the sum of the products is 

Ian"] X + [6n] y + [cn\ z + {dn} w -|- [nn"] = [vn] = [vv] 

If this equation be annexed as a fifth normal equation to the 
group (66), and the successive substitutions are made in it as in 
the others, beginning with x, it evidently becomes, successively, 

[6n.l]y -f [cn.l]2: + ldn,l]w + [nn.l] = [vr] 

[en . 2] 2r + Ldn . 2] m? + [wn. 2] = [tw] 

[dn.3]M? + [nn.3] = [vr] 

[nn.4] = [w] 

which last is the same as (71). 

DETERMINATION OP THE WEIGHTS OF THE UNKNOWN QUInTITIBS 
WHEN THE ELIMINATION HAS BEEN EFFECTED BY THE METHOD OF 
SUBSTITUTION. 

45. By the general method explained in Art. 86, the elimina- 
tion would have to be performed as many times as there are 
unknown quantities. It is desirable to have more direct methods. 
When there are but four unknown quantities, we can find their 
weights from the auxiliaries occurring in two successive elimina- 
tions in inverse order. In the first elimination, according to the 
order a, 6, c, dj we find w by substitution in the last normal 
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equation, and, the coefficient of to being then [dd . 3], it follows, 
by Art. 36, that the weight of the value of lo is 

p^ = ldd,S:\ 

In the inverse elimination, in the order d, c, 6, a, the coefficient 
of X in the final equation, which would be denoted by [aa . 3], 
will be the weight of x, or 

i>, = [««-3] 

Now, if a third elimination were carried out in the order x, y, w, z, 
or a, by rf, c (the third normal equation now taking the last place), 
we should have the same auxiliaries as in the first elimination, 
BO far as those denoted by the numerals 1 and 2 ; and the equa- 
tions (68) would still be the same, but in the following order : 

ldd,2'] w + lcd,2'] z + [rfn.2] = 
\cd.2'\w + [cc.2]2: + [cn.2] = 

The value of w given by the first of these is 

[0^.2] ^^ [(fn.2] 
\dd.2:\ [dd.2] 

which, substituted in the second, gives for the coefficient of z^ 

Therefore we have 

\dd . 3] 



i), = [cc.2] 



[(^.2] 



In the fourth supposed elimination, in the order d^ c, a, 6, the 
auxiliaries denoted by 1 and 2 would be the same as in our 
actually performed second elimination ; but in the final equation 
in y we should have for the coefficient of y the quantity 

[66.3] = [W.2] - M:|l [a6.2] = [aa.8] X ^ 
and, therefore, 

» =[66. 21 

Thus, when the elimination has been once inverted, we have 
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found the weights of two of the unknown quantities directly, 
and the weights of the other two in terms of the auxiliaries pre- 
viouslji used, and in a form adapted for logarithmic computation. 

46. In order to give the above method greater generality, so 
that the reader may be enabled to extend it to a greater number 
of unknown quantities, we remark that the product of the form 

P=laa] [66.1] [cc.2] Idd.S] 



has the same value whatever order may be followed in the elimi- 
nation. This is the ^ame as saying that it is a symmetrical func^ 
tion of a, 6, c, rf . . . which is, consequently, not affected in value 
by the permutation of these letters.* Suppose, then, four ordere 
of elimination, in which each unknown quantity in turn becomes 
the last, while the order of the remaining three quantities 
remains the same ; and, to distinguish the auxiliaries which occur 
in each elimination, let the letter which occurs in the last auxiliary 
be annexed to each of the others ; the above constant product 
may thus be expressed in the following four forms : 

P = lad]^ [66.1]^ [^^ .2]^ Idd.S] 

= [««]e [^^ • l]c [^ • 2]^ ICC . 3] 

= [aa\ Ice . l]j [dd . 2]^ [66 . 3] 
= [66]. [cc.l].[(i(i.2].[aa.3] 

Now, it is evident that each time a new unknown quantity is 
made the last, we do not change all the auxiliaries, but only 
those which involve the letter which has become the last in the 
new order. It is readily seen, therefore, that if we annex a letter 
to those auxiliaries only which have a different value from that 
which is denoted by the same symbol in the first elimination, we 
shall have, simply, 

P = [aa] [66 . 1] [cc . 2] [dd.^'] 
= [aa] [66.1] ldd,2] [cc .8] 
= lad] Ice .1] Idd. 2\ [66 . 3] 
==[66][cc.l].[(W.2].[aa.3] 

* The quantity P is, in fact, nothing more than the common denominator of the 
Talues of Xf y, z, w, when these values are reduced to functions of the known quan- 
tities and in the form of simple fractions ; and this common denominator must eyi- 
dently hare the same value whatever order of elimination is followed. 
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from which we deduce 



[rfd.3] 



^[6i.3] 



P. - [«' 



[c*;.l] [dd.2]. 



-- [aa] ■ 



jbb . 1] [CC.2] [rf d.3] 



If this method ia applied in the case of six unknown quantities, 
we shall in each of two eliminations have the weights of three 
of the unknown quantities by computing each time but one new 
auxiliary, and, therefore, the weights of all six when the second 
elimination is the inverse of the first. In the case of but four 
unknown quantities, by inverting the elimination we can find 
the weights of 2 and y twice, and thus verify our work. 

47, If we have but three unknown quantities, the weights are 
determined at the same time with x, y, and z themselves, by a 
single eliminatiou in the order a, b, c, in which z comes out first 
with the weight 

?. = [«:. 2] 

and then y and z, with the weights 
^^H j,^ = [bb . 2] = [bb . 1] 

1 J5. — \aa . 2] = [aa] 



in which 



[ec . 1], = [cc] 




INDEFESDENT DETERMINATIONj^OP EACH UNKNOWN QUANTITY AND 
ITS WEIGHT, ACCORDING TO SADSS. 

48. Let the four equations (70) he multiplied respectively 
by 1, A', A", A'", and let these factors be determined by the 
condition that in the sura of the products the coefficients of y, 
z, and w shall be zero. Also, let the last three equations of (70) 
he multiplied respectively by 1, B", B'", aud let these factors 
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be determined by the condition that in the sum of the products 
the coefficients of z and w shall be zero. Finally, let the last 
two equations of (70) be multiplied respectively by 1, C'"^ and 
let C"' be determined by the condition that in the sum of the 
products the coefficient of w shall be zero. The conditions 
which determine these factors are then 

= \^ + A' 

[aa] ^ 

[aa] ^[66.1] ^ 

[aa] "^[66.1] "^[cc. 2] "^ . ^^^^ 

[66.1] ^ 

[66.1] "^[cc. 2] "*" 

[cc.2] "^ 

and the final values of x, y, z^ w^ in terms of these factors, are 
given as follows : 

_ ar = ^ + t^^^ ^ ^' 4- [g^'2] ^f/ . [^^-^] j^/// 
[aa] "^[66.1] "•" [cc.2] "•" [tW.S] 



^ [6n.l] [en. 2]^,, [^n . 3] ^,,, 
^ [66 . 1] "•" [cc . 2] "^[(Zd.3] 



[cc.2] "^ [(id. 3] 
[(Zn . 3] 



(74) 



[dd,'^'] 

49. As the equations (73) are above arranged, all the factors 
A are determined from the first system of three equations ; the 
factors B from the second system of two equations, &c, ; in each 
case, by successive substitution. This method then enables us 
to find each unknown quantity independently of the others. 

Another form may be given to the computation of the auxiliary 
factors. Since in the formation of the equations (74) we have 
regarded [an], [6/i], [en], &c. as independent, we must still so 
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them when we invert the process and recompose the 
ions (70) from (74). I^ then, we moltiply the equations (74) 

jtively by 1, p— i ^_J i_J, and add the products in order 

lad] laa] [ad] ^ ^ 

recompose the first of (70), the coefficient of [anl will be - — j^ 

it the coefficients of [bn . 1], [m . 2], &c. must severally be equaJ 
zero. The same principle will apply when we recompose the 
md equation of (70) from the last three of (74), &c Hence 
Fwehave 

[oft] 



= ^' + 



[oa] 
[oft] 



« = ^"+H^'' + 



{ae2 
lad] 



0=.A'"+f^B'" + 



[^*] T>/// , [^]gw 



= J?" + 



0=B'" + 



= C" + 



laa] 

Ibc.l] 

[W.l] 

Ibc.l] 

Ibb.l] 

led. 2] 

Ice. 2] 



laa] 



+ 



lad] 
laa] 



(75) 



C'" + 



Ibd.l] 
Ibb.l] 



According to this scheme, we first find A\ -B", C" from the 
equations in which they occur singly ; then, with these factors, 
we find the values of -4", -B'", from the equations involving two 
&ctors, &c. 



50. Again, let us write the 3d, 5th, and 6th equations of (75) 
in the following order : 



^"/+[^^^'^+M(7'" 



laa] 



B'" + 



laa] 

Ibc.l] 

Ibb.l] 



+ 



C"" + 



C'" + 



[ad] _ 
laa] 
Ibd.l] 
Ibb . 1] 

led. 2] 
Ice . 2] 



= 



= 



= 



Comparing these with the first three of (70), we at once infer 
that J.'", -B'", C" are those values of x, y, Zj respectively, which 
we should obtain from our first three normal equations by putting 
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w =-1 and omitting the terms in n ; or, going back to (66), tliat 
-4'", -B'", C"' may be determined by the following conditions: 

[aa] A'" + [oft] B"' + [oc] C" +[ad] = 
[aft] A"' + [66] J5'" + [6c] C" + [6d] = 
lac^ A''' + [6(?] -B'" + [cc] 0"' + [cd] = 

If now we multiply the normal equations (66) by J.'", -B'", C", 
and 1, respectively, and add the products, the conditions just 
given will cause Xy y, and z to disappear, and the resulting equa- 
tion in w must be identical* with (69) : so that A''\ B'", C" 
must also satisfy the following condition : 

[an] A'" + [6n] J5'" + [cn\ C"' + [dn\ = [dn . 3] (76) 

The second and fourth equations of (76) being written as follows, 

\aa\ [aa] 

[^cj]_ 
^ [66.1] 

and compared with the first two of (70), we infer that -4", 5" 
are those values of x and y which we obtain from the first two 
normal equations by putting z = 1, 2(7 = 0, and omitting the 
terms in n; that is, A" and B" must satisfy the conditions 

[aa] A" + [a6] B" + [ac] = 
[a6] A" + [66] B" + [6c] = 

Therefore, if we multiply the first three normal equations (66) by 
A'\ JB'', 1, respectively, and add the products, x and y will dis- 
appear, and, the resulting equation being identical with the first 
of (68), we must also have 

[an] A" + [6n] 5" + [crq = [en . 2] (77) 

Lastly, it is evident that A' must also satisfy the condition 

Ian] A' + [6n] = [6n . 1] (78) 

From these relations we readily infer general formulae for the 

weights of the unknown quantities. 

. — I 

* The equation (69) is the last normal equation, unchanged except by tlie substitii- i 
tion of equivalents for z, y, and z ; and in the present article we eliminate z, jr, and t 
by the use of factors, but do not change the last normal equation, since we multiply 
it by unity. 
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' According to Art 34, the reciprocal of the weight of z is that 
value which we obtain for x if we put [an] == — 1 and [6n] = [en] 
= [c/w] = 0. But, under these conditions, the equations (76), 
(77), (78) give 

[e^n.3] = — A'", [en. 2] = — -4", [6n.l] = ^A' 

In order, therefore, that the value of x given by the first equa- 
tion of (74) may become — , we have only to substitute — -4'", 

— -4", — J.', — 1, respectively, for [dn , 3], [en . 2], \bn . 1], [an]. 
In the same manner, the weight of y being found by putting 
[6n] = — 1 and [<m\ = [en] = [dn] = 0, we have to put 

[(in.3] = — ^", [cn.2] = — ^', [6».1] = — 1 

in the second equation of (74), in order that we may put — for y. 
For the weight of z we have to put 

[dn.3] = — C", [en. 2] = — 1 

in the third equation of (74), and — for z. 

I'm 

For the weight of -w?, we have to put 

[(Zn.3] = — 1 
in the last equation of (74), and change w to — • 
The final formulae for the weights are, therefore, 

11. A'A' , A" A" . A''' A''' 



i>, laa] [66.1] [ce.2] ' [dd.S} 

1 1 B"B" W'B'* 



(79) 



i>^ [66.1] ' [ee.2] \M.%-\ 

l>,'~[ee.2]"^ [<W.3] 
1 1 

MEAN ERROR OF A LINEAR FUNCTION OF THE QUANTITIES X, y, 2, W. 

50. To find the mean error of ttie function 

X = fx + jfy + A^ + IM? + / (80) 

xohjm Xy y, Zj w are dependent upon the same observations. 



542 APPENDIX. 

The quantities a:, y, z^ w not being directly observed, their 
mean errors cannot be treated as independent, as was done in 
the case of directly observed quantities in Art. 22. "We might 
proceed by the method of Art. 28; but, as we here suppose 
X, y^ Zy w to have been determined from the normal equations 
(66), we can obtain a more convenient method by the aid of the 
auxiliaries which have been introduced in the general elimina- 
tion. The quantities x, y, ^, w being functions of the directly 
observed quantities n', n", w'", . . . the mean error of J^ can be 
readily obtained by the principles of Art. 22, if we first reduce 
X to a function of these observed quantities. For this purpose, 
if the values of x, y, z, w deduced from (70) be substituted in X^ 
we shall have an expression of the form 

X= k^ [an] + k^ [6n.l] + A, [en. 2] + A, [rfn. 3] + I (81) 

in which the coefficients k^, Aj, k^j k^ are functions of [aa], [aft], 
&c. In order to determine these coefficients, let us substitute in 
this expression the values of [an], [6n.l], &c. given by (70). We 
find 

Jr= — lad]kf;c — [^l^oV — [^]^o^ — [a^Afo^+i 

— U>b .l]kjf --[be A]k^z -^[bd .l] k^w 

— [cc . 2] Z:^ — [cd . 2] k^w 
— [<W . 3] /:,u? 

which becomes identical with (80) by assuming 

Jiad]k^=—f 

[a6]A:,+ [66.1]/:, = ~(7 

lac^k^ + [bcA]k^+lcc.2]k^ = — h ( C»^> 

[a^q k^ + Ibd.l] k^ + [c(Z.2] k^ + [dd.3] /:, = — i 

These equations fully determine the coefficients. We find k^ 
directly from the first, and then A^, k^y Aj, by successive substitu- 
tions in the others. 

Now, to find the mean error of JT under the form (81), let the 
mean error of each of the observed quantities w', w", n"' .... be 
denoted by e (these observed quantities being supposed of eqniJ 
weight, or, rather, the equations of condition being supposed to 
have been reduced to the same weight), and let the correspond- 
ing mean errors of 

lan]f [ftn.l], [<?n.2], [rfn.3], -T, 
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be denoted by 

S. E 
Siuce we bave 


-^« 


B., (.JO. 


[a»] = «■ 


»■+«"»■ 


+ a"V"+.... 


we have, by Ajrt 
Again, we bave 

and henec 


22, 

= M - 






K 


^([46] 




W + glM) 


P 


h([»6] 




= 


[6s.i:^ "■"■■ 




Iq a Bimilar man 


ner, we 


lave, also. 


-«. 


= [cc. 


2].', 


i:,'=ldd.31,< 



The quantities x, ?/, 2, ?f, being determined from the equations 
(70), their mean errore involve those of the quantities [an], [bn . 1], 
[en. 2], [dn.S], precisely as if the latter had been independently 
observed quantities affected by the mean errors just determined. 
Hence also in (81) we regard [tm], [6n.l], &e. as independent; 
and it then follows directly from the principles of Art. 22 that 



(eJO' 



= (V [««] + -^i' l'>'> ■ 1] + V ["" ■ 2] + ft.' ['i'^ ■ 3]) E^ 



51. Prom the preceding article we may easily find the for- 
mula (74) and (79). The function ^ becomes x when we assume 
/^^Ij? — ^ — *'=^ — 0; and then (81) gives x while (83) gives 
e^', and hence tiie weight ^ — ,■ This hypothesis gives in (82) 
[ad] A'u = — 1 ; and the remaining equations of (82) are identical 
with the first three of (78) if we put [66.1] k, = —A', [c(t.2] i, 
= — A", Idd. 3]ij = — A'"; andthen (81) becomes identical with 
the first of (74), and (83) with the first of (79). la a similar 
manner we may deduce the remaining equations of (74) and (79). 
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Example. — In order to exliibit the numerical operations which 
the preceding method requires, in their proper order and within 
the limits of the page, I select an example involving but three 
unknown quantities. The following equations of condition were 
proposed by Gauss {Theoria Motus Corp, Coel.y Art. 184) to illus- 
trate his method: 

(1) X— y + 2z= S 

(2) 8a: + 2y — 5^= 5 

(3) 4a: + y + 4^ = 21. 

(4) — 2a: + 6y + 62: = 28 

of which the first three are supposed to have the weight unity, 
while the last has the weight J. Multiplying the last by |/ J = J 
(Art. 41), the equations of condition, reduced to the same weight, 

are — 

(1) a; — y -f. 22: — 3 = 

(2) 3a: + 2y — 5-? — 5 = 

(3) 4a: + y + 4^: — 21 = 

(4) — a: + 3y + 32: — 14 = 

The next step is to form the coeflGLcients [aa], [aft], &c., of the 
normal equations. In the present example this can be done very 
easily without the aid of logarithms ; but, in order to exhibit the 
work usually required in practice, I shall give the forms for 
logarithmic computation. The sums of the coefiSLcients of the 
unknown quantities will be employed as checks, according to 
Art. 80. Their logarithms, together with those of a^bjC^iiy are 
given in the following table : 



log a 


log b 


logc 


log« 


logn 


0.00000 


nO.OOOOO 


0.30103 


0.30103 


n0.47712 


0.47712 


0.30103 


nO.69897 


00 


nO.69897 


0.60206 


0.00000 


0.60206 


0.95424 


nl.82222 


fiO.OOOOO 


0.47712 


0.47712 


0.69897 


nl.l4618 



(1) 

(2) 
(3) 
(4) 

It is important, where many operations are to be performed, to 
write down no more figures than are necessary for the clear prose- 
cution of the work. Hence, in combining the preceding 
logarithms it will be found expedient to proceed as follows. 
Write each log a upon the lower edge of a slip of paper; then, 
placing this slip so that log a shall stand over log a, log 6, log e, 
&c., of the same horizontal line, in succession, add together tibe 
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two logarithms jnentaliy, and, with the Bum in the head, take from 
tlie logarithmic table the corresponding Datural number (aa, ab, 
ac, as, or an], which place iu a column appropriated for the pur- 
pose. Theu write log b in the same manner, and form bb, be, bs, 
bn, and so proceed to form all the coefficients of the normal 
equations, as iu the following table : 





[aa] 


l-"} 


[" 


=J 




['' 


'1 


[a 


n] 


m 


t^] 1 




+ 


+ 


— 


+ 




— 


+ 


— 


+ 


^ 


+ 


+ 


— 






















R( 


l.C 




2.0 


"1 


9.0 


fi.f 








Ifil 




().(: 




Ifil 






10.0 


h) 


16.0 


41 




Iflt 






am 






B4f 


1,( 


4.C 




(*) 


1.0 




a.o 






JLU 




5.0 


14.0 




9.0 


9.0 








ll».t 


41 


IKI 




IMI 


88.01 G-t 


14.1 


Wi.\ 










+ 27.0 


+ 


.u 










+ 


iJ.U 




+ 1G.0 


+ 







[h] 


[bn] 


[«] 


["1 


i-: 




[«J 


f-l 




+ - 


+ 


— 


+ 


+ 


— 


+ 




+ 


— 


+ 


1) 


2.C 


sr 




4.r 


4.1 






fir 




6,i: 


9.0 


0.( 




[IH 


25.C 




fl.r 


•/hi 




0,( 




2&.0 








■'1 1 


18.C 


KHI 






K4I 




TRii.r 


441.0 


41 


15.( 




42.0 


9.0 


16.U 






42.1 




70,0 


196.0 








v:ii 




W. 1 


IK. 


■:'h ( 


\H-i.l 










+ 22.0 




0.0 


+ 54.0 


+ 


5.(1 


— 1 


07.0 


— 2 


66.0 


+ 671.0 



Having ascei-taincd that the results eatiBty the test equations 
(48), we can write out the normal equations as follows: 



7x 


+ 


6j 


— 


88 


lix 


+ 


16v + 


s — 


70 






» + 


54s — 


107 



= , 

"We proceed to determine the values of x, y, z, according to 
our general formulse, still carrying out the work with logai'itbms 
foi' the sake of illustration. Here, again, sjstem and concise- 
ness are indispensable. The whole computation is given below 
nearly in the form proposed by Encke. This form corresponds 
to the group of equations (70). It is divided into three principal 
compartments, corresponding, respectively, to the first three equa- 
tions of (70), each beginning one column farther to the right. In 
the first compartment the first line of numbers contains the values 
of [(m], [«i], kc, the second line their logarithms, and the third 
line the logarithms of the coefficients of the first equation. The 
logarithms in this third line are formed by subtracting the first 
log, in the second line from each of the subaequent ones, for this 
Vou 11.-35 
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purpose writing the first logarithm upon the lower edge of a slip 
of paper. 

In the second compartment, the first line contains the values 
of [bb], [6c], &c. ; the second line, the quantities subtractive from 
these, according to the formulae in Art. 42. To form these sub- 
tractive quantities, write the logarithm of ^ — ^ (which is here 

9.34679) upon the lower edge of a slip of paper, and hold it suc- 
cessively over log [ah"] and each of the subsequent logarithms in 
the same line ; add the two logarithms mentally in each case, take 
the corresponding natural number from the logarithmic table, 
and write it in its place below. Subtracting these numbers, we 
have the values of [W.l], [6<?.l], &c. The fourth line contains 
the logarithms of these quantities ; the fifth, the logarithms of the 
coefficients of our second equation, formed by subtracting the 
first logarithm of the preceding line from each of the subsequent 
ones in that line. 

In the third compartment we have — first, the values of [cc], &c. ; 
secondly, the values of the subtractive quantities formed from 
the last line of the first compartment as before ; thirdly, the 
remainders which are the values of [cc.l], &c. The fourth line 
contains the values of the quantities which are subtractive from 
the preceding and are formed from the last line of the second 
compartment by adding the first logarithm of that line to the 
logarithm immediately above it and to each of the subsequent 
logarithms in the same line ; the fifth line contains the remain- 
ders which are the values of [cc.2], &c. ; the sixth line, the loga- 
rithms of these ; and the last line, the logarithms of the coeffi- 
cients of our third equation. 

For control, we carry through the operations upon [as]^ [b$], 
&c., precisely as upon the other quantities; and then, according 
to the arrangement of the scheme, we should have, if we have 
computed correctly, each sum containing s equal to the sum of 
the quantities on its left in the same line, together with those of 
the same order in a vertical column over the first number in this 
line. Thus, we must have, in the present case, 

lbs. I] = [bb . 1] + [6c. 1] Isn . 1] = [6n . 1] + [en. 1] 

[cs.l] = [cc.l] + [6c.l] [5n.2] = [cn.2] 

[cs . 2] = [cc . 2] 

relations easily proved by means of the formulae of Art. 42 com- 
bined with (48). 
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The columns [sn] and [vri] are added to the third compart- 
ment in order to form the quantity [wn . 3], from which the mean 
error of observation is to be deduced, as will be shown hereafter. 



[aa] 

+ 27.000 
1.43136 



[ab] 

+ 6.000 
0.77815 
9.34679 



[ac] 
0.000 

— QO 
00 



— 88.000 

0.000 

-h 21.306 



— 66.695 
nl.82409 
logx = 0.39273 



[bb] 

-f 16.000 
+ 1.333 



-f- 13.667 
1.13666 



— 50.444 

4- 1.916 

— 4^528 
nl.68599 

logy=: 0.55083 



[be] 

+ 1.000 
0.000 



-f- 1.000 
0.00000 
8.86434 



[cc] 

-f- 64.000 
0.000 



-h 54.000 
4- 0.073 



+ 53.927 
1.73181 



[CM] 

+ 83.000 

1.51851 
0.08715 



[bs] 

+ 22.000 
+ 7.333 



+ 14.667 
1.16633 
0.03067 



+ 56.000 
0.000 



+ 55.000 
+ 1.073 



+ 63.927 



[an] 

— 88.000 
nl. 94448 
ftO.51312 



[6n] 

70.000 
19.556 



- 50.444 
111. 70281 
wO.56715 



[en] 

— 107.000 
0.000 



— 107.000 
- 3.691 



— 108.309 
n2.01414 



[m] 

— 265.000 

— 107.555 



— 157.445 

— 54.135 



— H>8.810 



[nn] 

+ 671.000 
+J86.813 

+■384.187 
+ 186.191 

+197^0 
+ 197.901) 



log (_ z) = nO.28283 [nn . 3] = + 0.087 



After z has been found, its value is substituted in the second 
equation of (70), and y is deduced. Then, the values ofy and z 
being substituted in the first equation, we find x. The numerical 
computations are given above in the margin. 

Then, for the weights, by Art. 47, we have first to find the 
additional auxiliary 

[cc.l].= [cc]-^][6c] 
and by the formulse of that article we have — 



m 


m 


+ 15.000 


+ 1.000 


5 1.17609 


0.00000 




8.82391 




Ice] 




+ 54.000 




+ 0.067 


[ce.la] = 


: + 53.933 



log[W.l] 1.13666 
log [66] 1.17609 


log lee . 2] 
log [cc . 1] 
log [cc . 1] a 


1.73181 

1.73239 

I 1.73185 


1.43136 
9.95957 
9.99996 


1.13566 
9.99942 


1.73181 
logl>. 


1.13508 




1.39089 
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The final iresult is then 

x=: + 2.4702 with the weight 24.597 
y = + 3.5508 ** « 13.648 

z = + 1.9157 « « 63.927 

It only remains to substitute the values of Xy y, and z in the 
original equations of condition, to form the residuals v, and from 
these to determine the mean error of observation. Since here 
there are but three unknown quantities, we have, by (71), 

[vv] = [nn . 3] 

and hence the mean error of an observation of the weight unity 
is, by (61), m being the number of equations of condition, 



=i'^>)='-^ 



The direct computation of the residuals is, therefore, not necessary 
for determining e: nevertheless, it is desirable in most cases to 
resort to the direct substitution also, not only for a final verifica- 
tion, but in order to examine the several observations, and to 
obtain the data for rejecting any doubtful one by the use of 
Pbirce's Criterion, to be given hereafter. This direct substitu- 
tion has already been carried out for this example on p. 525, 
where we have found [vv] = 0.0804, which agrees with the above 
value of [nn. 8] as nearly as can be expected with the use of five- 
decimal logarithms. 

52. It not unfpequently happens that one of the unknown 
quantities is such that the given observations cannot determine 
it with accuracy. For example, in the reduction of a number 
of observations of an eclipse, one of the unknown quantities is a 
correction of the moon's parallax ; but, unless the places of ob- 
servation be remote from each other, the correction will be very 
uncertain, and this uncertainty will affect all the other quantities 
which enter into the equations of condition. In such a case, this 
unknown quantity will come out with a small coefficient, which 
of itself will reveal the existence of the uncertainty when it is 
not otherwise anticipated. In order that this uncertainty may 
not affect those quantities which are well defined by the obser- 
vations, it is expedient to determine all the latter as functions of 
the uncertain quantity, which for that purpose must be made the 
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last in the elimination. Thus, "W'ith four unknown quantities 
T, y, 2, ic, we proceed only as far as the auxiliariee denoted by 
the numeral 2; then, having found the factors A', A", A'", j5". 



I 



,bj(73)t 



(75), if we put 




[ST] + [SITT] ^ + [< 


t.2J 


[M.l] 1 L«.3J 




[«.2J 





these will give the values of the unknown quantities which we 
should obtain from the first three normal equations if the last 
unknown quantity were disregarded or pat = 0. Then, by (74), 
the final vahies of x, y, z, as functions of the uncertain quantitj' 
w, will be 

a: =^ ic' -f A"'w ^ 

y = y- + Sf"v> \ (85) 

z =if -\- C"'w ) 

The values o(x', y', ?', will thus be well determined, and a sub- 
eequent independent determination of w will enable us to find 
the final values of a;, y, z.* 

Having found the weights of x', y', z' (which is done as if they 
were the only quantities under consideration), and their mean 
errors e^'. e^', e^', then, when the quantity w is afterwards fouiid, 
the mean errors of the final values will be 

■■.= ■.■■+( J"'0' ■ 

■*,= ■,'•+ (-B"'0' 

•*.=■."+(<;"'..)■ 

as we find from the equations (79), or by Art. 20. 



\ qTia: 



COSDITIONED OBSERVATIONS. 
18. In all that precedes, we have supposed that the 



lantities to be found by observation, either directly or indirectly. 
were independent of each other. Although they were required 
to satisfy certain equatioufi of condition as nearly as possible, yet 
they were so far independent that no contradiction was involved 
in supposing the valuesof one or more ot'them to be varied without 
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varying the others. By such variations we should obtain sys- 
tems of vahies more or less probable, but all possible. 

There is a second class of problems, in which, besides the 
equations of condition which the unknown quantities are to 
satisfy approximately, there are also equations of condition which 
they must satisfy exactly: so that of all the systems of values 
which may be selected as approximately satisfying the first kind 
of equations, only those can be admitted as possible which satisfy 
exactly the equations of the second kind. The number of these 
rigorous equations of condition must be less than the number of 
unknown quantities ; otherwise they would determine these 
quantities independently of all observations. These rigorous 
equations, then, may be satisfied by various possible systems of 
values, and we can therefore express the problem here to be con- 
sidered as follows : Of all the possible systems of values ivhich exactly 
satisfy the rigorous equations of condition, to find the most probable, or 
that system which best satisfies the approximate equations of condition. 

The following are simple examples of conditioned observations. 
The sum of the three angles of a plane triangle must be 180°: so 
that if we observe each angle directly, and the sum of the observed 
values differs from 180°, these values must be corrected so as to 
satisfy this condition. The sum of the angles of a spherical 
triangle must be 180° + spherical excess. The sum of all the 
angles around a point, or the sum of all the differences of azimuth 
observed at a station upon a round of objects in the horizon, must 
be 360?. 

The approximate conditions in these cases are expressed by 
the observations themselves ; for the final values adopted must 
correspond as nearly as possible to the observed values. The 
corrections to be applied to the observed values are to be re- 
garded as residual errors with their signs changed ; and the solu- 
tion of our problem is involved in the following statement: Of 
all the systems of corrections which satisfy the rigorous equations, that 
system is to be received a^s the most probable in which the sum of the 
squares of the residuals in the approximate equations is a minimum. 

54. The general problem as above stated may be reduced to 
that of unconditioned observations, already considered. For let 
us suppose there are m' rigorous equations of condition, and m 
unknown quantities. From these m' equations let the values of 
m' unknown quantities be obtiiiiied in terms of the remaining 



METHOD OF LEAST SQUARES. 



551 



m — m' quantities, and let these values be substituted in all the 
approximate equations of condition; then there will be left in the 
latter only m — m' quantities, which may be treated as independ- 
ent, so that, the approximate equations being now solved by the 
method of least squares, we have the values of the m — m' quan- 
tities, with which we then find the values of the first m' quan- 
tities. This is a general solution of the problem; but it is not 
always the simplest in practice. I shall illustrate it by a simple 
example, before giving a method applicable to more complicated 
cases. 

Example. — ^At Pine Mount, a station of the U. S. Coast Survey, 
the angles between the surrounding stations 1, 2j 3, 4 were 
observed as follows : 



1.2 
2.3 
3.4 
4.1 



Joscelyne— Deepwater 65^ 11' 52".500 

Deepwater— Beaky ne 66 24 15 .553 

Deakyne— Burden 87 2 24.703 

Burden— JoBcelyne 141 21 21 .757 



weight 
3 
3 
3 
1 



There are here four unknown quantities subjected to the single 
rigorous condition that their sum must be 360°. But, instead of 
taking the angles themselves as the unknown quantities, we shall 
assume approximate values of them, and regard the corrections 
which they require as the unknown quantities. 
We assume 



1.2 
2.3 
3.4 
4.1 



Joscelyne — Deepwater, 65° 11' 52".5 + w 
Deepwater — Deakyne, 66 24 15 .5 + ^ 
Deakyne—Burden, 87 2 24 .7 + y 

Burden — Joscelyne, 141 21 21 .8 + 2 



the sum Of which must satisfy the condition 

359° 69' 54".5. + u? + a: + y + 2 = 360° 
or 

w -\-x + y -{-z^ 5".5 == 

The difference between the assumed value and the observed 
value in each case gives us a residual; and the approximate 
equations of condition are, tiierefore^ 

w— =0 

X — 0.053 = 

y _ 0.003 == 

' z -\- 0.043 = 



552 



APPENDIX. 



We have here but one rigorous condition (or m' = 1), and to 
eliminate this we have only to find from it the value of one un- 
known quantity in terms of the others, and substitute it in the 
approximate equations of condition: thus, substituting the value 

w = — X — y — 2 + 5".5 

our equations of condition, containing now three independent 
unknown quantities, are 

weight. 



— a: — y — 2 + 5".5 = 

X — .053 = 

y —0 .003 = 

2: — .043 == 



3 
3 

1 



The normal equations, applying the weights, are then 

6x + 3y + 3^ — 16.659 = 
3x + 6y + 3^ — 16.509 = 
3a; -^ 3y + 4^ — 16.457 = 



which, being solved, give 



whence also 



x = + 0".9675 
y = + .9175 
2 = + 2 .7005 

v> = + 0.9145 



and the corrected values of the angles are 



1.2 
2.3 
3.4 
4.1 



Joscely no— Deep water 65** 11' 53".4145 

Deep water— Deakyne 66 24 16 .4675 

Deakyne— Burden 87 2 25 .6175 

Burden— Joscelyne 141 21 24 .5005 

360 .0000 



55. "When the number of unknown quantities is great, or when 
there are several rigorous conditions to be satisfied, the preceding 
method would lead to very tedious computations, since we are 
required to perform two eliminations, the first from our m' 
rigorous equations to find the first m' quantities in terms of the 
others, and the second from our normal equations involving all 
the remaining quantities. In order to obtain the general form 
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for a more condensed process, let the most probable values of a 
number (m) of directly observed quantities be 

F', F", F'", &c. . . . F<«' 
Let the observed values be 

M', M'\ M"', &c. . . . Jlf <•> 
Let these observations have the weights 

jp', fj f\ &c. ...jp<«> 

Let the equations which the most probable values are required 
to satisfy rigorously be expressed by 

^' =/' (F',F",F'",...) = 

^n ^ fii ( yi^ yff^ yw^ . . .) = 

^m ^ fin ( yi^ yn^ ym^ . . .) = ( C^"^) 

&C. 

and let 

m! = the Dumber of these conditions. 

Let the most probable corrections of the observed values be 

t/, t/', t/", &c. . . . !;<•> 
so that 

F' = Jf' + i/, F" = Jf " + t/', F'" = Jf' " + v'", &c. 

Let the values of ^', ^", f>'" . . . when the observed values are 
actually substituted be n\ n", w"' ... or 

/' (M', M'\ M'", ...) = n' 
/" (M'y M'\ Jtf"% . . .) = n" 
/'" (Jf' , Jf ", Jf' ", . . .) = n'" ' ( (^^) 

&c. 

Le. «.e diiTe^BtW oo,<Eci,u„ ^. ,^, *o., |^. ^L *o. be 

formed ; substitute in them the values M'^ M^'y M''' ... for F', 
F", F''', and denote the resulting values by a', a'', &c., 6', 6", 
&c. ; that is, put 

ifl = a', -^ = a", -^ = a'", &c. 
dF' dF" dF'" 

^ = 6', -^?!- = 6", ^?!- = 6"',&c. 
dV ' dV" dV" ' 

dV ' dV" ' dV" 
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These values of the differential coefficients will generally be suf- 
ficiently exact ; but if J!f' , ilf ", M*" ... are found very greatly in 
error, a repetition of the computation might be necessary, in 
which the more exact values found by the first computation 
would be used. 

The values of J!f' , M"^ M'" . . . being assumed to be so nearly 
correct that the second and higher powers of the corrections i?', 
v", v'" . . . may be neglected, we have at once, by Taylor's 
Theorem, as in the similar case of Art. 40, 

^' = n' + aV + a'V + a"V"+ . . . . + a^-V> = 

&c. &c. 

which m' equations must be rigorously satisfied by the values of 

The equations 

F' — Jf' = 0, F" — M" = 0, F'" — Jf' " = 0, &c. 

are the approximate equations of condition ; or, more strictly, 

F' — -af' = t;', F" — Jf"=t/', F" — -af '" = v"', &c. 

are the equations of condition which are to be satisfied by the 
most probable system of residuals t?', i?", v'" . . . .' These, reduced 
to the unit of weight by Art. 41, become 

( F' — M') i/p' = vyp^, ( V" — M") |//' = v'Yp'', &c. (90) 

and the most probable residuals v'|/p', v"\/p" are those the sum 
of whose squares is a minimum, or we must have 

yi;'* + y I?"' + />'"v"" + &c. = a minimum. 

Putting, then, the differential of this quantity equal to zero, we 
have 

p'vfd'if + p'^'i^'dv" + p'"v"'dv'" + &c. = (91) 

If v', I?'', v'" . . . were independent of. each other, each coeffi- 
cient of this equation would necessarily be zero (as in Art. 28), 
and then the most probable values of F', F", V" . . . would be 
the directly observed values M'^ M*\ M'" . . . But this minimum 
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ia here conditioned by the equations (89). If, then, we differen 
tiute (89), the equatious 

a'dv' -\- a"iW -f a"'dv"' -|- =; \ 

b'do' + h"dv" + y'dif" -f . . . . = / 

ddv' + c"di}" + if"dv"' + .... = ) (^^) 

&c. I 

must coexist with (91). 

The uumber of the equations (92) is m', while the nnmher of 
differentiaia is m; and since, by tlie nature of the case, we must 
have m > m', wc can, by elimination, find from (92) the values 
of m' differentials in terms of the remaining m — m' differentials. 
Let ua suppoae this elimination to be performed, and that the 
values of the first m' differentials, found in terms of the others, 
are tlieu substituted in (91) ; we ahall thus have an equation in 
which the remaining m — m/ unknow^n quantities can be regarded 
as independent, and the coefficienta of these m — m' quantities 
in this final equation will then severally be equal to zero. We 
can arrive directly at the reault of such an elimination and sub- 
stitution as follows. Multiply the fii-st equation of (92) by A, the 
second by B, the third by C, &e., and also the equation (91) by 
— 1, and form the sum of all these products. Then, if A, B, 
C . . . are determined ao that m' differentials shall disappear 
from the sum (and they can be so determined, since it only 
requires m' conditions to determine m' quantities), the final 
equation obtained will contain only the m. — m' remaining differ- 
entials. But, the latter being independent, their coefficienta must 
also he severally equal to zero; and hence we have, iu all, the 
following m conditional equations : 



a' A + 6'B -\- dC +....— pV 

a"A -1- b"B -I- <f'C -\- p"xf 

a"' A -f b"'B + d"C H — y V 

&c. &c. 



= 
= 



(93) 



If we multiply the first of these by — > the aeeond by — ■ &c,, and 
add the products, we have, by comparison with the first equation 
of (89), 



[-].<+[?]-[-]« 
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in which the usual notation for Bums is followed. In this way 
we can form wl normal equations containing m' quantities, 
namely, 

[y]'' + [f]* + [f]<'+-+«'=» 

[f]^+[f]^+m''+-+'^=» ( w 



[f>+[T]^ + [y]<' + 



• • • • 



+ n'"= 



&;c. 



K the observations are of equal weight, we have only to put 
p = 1, or, in other words, omit p. 

The factors J., jB, C . . . are called by Gauss the correlatives of 
the equations of condition. 

The equations (94) being resolved by the usual method of 
elimination (Art. 42), the values of the correlatives found are 
then to be substituted in (93), whence we obtain directly the 
required corrections, 

1/ = -V («'^ +1>'B -^cfG +....) 

t/' = -1^ {a!' A + U'B + c"(7 + . . . .) , ,^,^ 

i/" = -4? («'"^ + V'^B + c^'C + . . . .) 
&c. &c. 

and hence, finally, the most probable values of the observed 
quantities, V = M' + v', V = JHf " + v'\ &c. 

The comparative simplicity of this process will best be shown 
by applying it to the example of the preceding article. We 
there have given, by observation, 

W = 65^ ir 52".500, / = 3 

M" = 66 24 15 .553, f = 3 

Jlf' " =87 2 24 .703, p"' = 3 

JIT^ = 141 21 21 .757, p'"" = 1 

with the condition 

F' + F" + r"' + P^ — 360^ = 
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We have, first. 

a' =: a" = a'" = a" = 1 

and when M', M", &c are put for V, V", &c., we have (88) 

»' = — 6".487 

Ab we have hut one condition, we have alao hut one correlative 
A ; the equation of condition is, by (89), 

— 5".487 + v' + v" + if" + i>'' = 

and the single normal equation may be conBtruct«d according to 
the following form : 



F 




p 


3 


1 


i 


3 


1 


» 


3 


1 


i 


1 


1 


1 


T 


h 


-2 



and hence, by (95), 

t/ = + 0.9145 
«" = + 0.9145 
v"'= + 0.9145 
v'' =4-2.7435 



2i — 5".487=0 
A = -i- 2".74S5 



Comctcd TftloH. 
V = 66" 11' 58".4145 
F" = 66 24 16 .4676 
F"'= 87 2 25 .6176 
F"= 141 21 24 .5006' 
S60 



agreeing -with the result found by the much longer process of 
the preceding article. 

56. The further prosecution of this branch of the subject 
belongs more especially to works on Geodesy. For more ex- 
tended examples, see the special report of Mr. C. A. Schott in 
the Report of the Superintendent of the U. 8. Coast Survey for 
1854, from which the above example has been drawn. Consult 
also Bbssel's Gradmeasang in. Ostpreussai m 1838; Rosehberqer, 
in the Astronomische Nachnchten, Nos. 121 and 122 ; Bessel, ibid. 
No. 438 ; T. Galloway, Application of the Method to a Portion 
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of the Survey of England, in the Memoirs of the Royal Astronomi' 
cal Society^ Vol. XV. ; J. J. BiEYBR's Kustenvermessung ; Fischer's 
Geodmsk; Gerling's AusgUkhungs JRechnungen ; Diengbr's J[w5- 
gUichung der Beobachiungsf elder ; Liagre, Galcul des ProbabMis; 
and Gauss, Supplemenium thecri(B combinaiionis, kc, 

CRITERION FOR THE REJECTION OP DOUBTFUL OBSERVATIONS. 

57. It has been already remarked (p. 490) that the number of 
large errors occurring in practice usually exceeds that given by 
theory, and that this discrepancy, instead of invalidating the 
theory of purely " accidental" errors, rather indicates a source 
or sources of error of an abnormal character, and calls for a 
criterion by which such abnormal observations may be excluded. 
The criterion proposed by Prof. Peirce* will be given here with 
the investigation nearly in the words of its author, and with only 
some slight changes of notation. 

58. "In almost every true series of observations, some are 
found which differ so much from the others as to indicate some 
abnormal source of error not contemplated in the theoretical 
discussions, and the introduction of which into the investigations 
can only serve, in the present state of science, to pei-plex and 
mislead the inquirer.^ Geometers have, therefore, been in the 
habit of rejecting those observations which appeared to them 
liable to unusual defects, although no exact criterion has been 
proposed to test and authorize such a procedure, and this delicate 
subject has been left to the arbitrary discretion of individual 
computers. The object of the present investigation is to produce 
an exact rule for the rejection of observations, which shall be 
legitimately derived from the principles of the Calculus of Pro- 
babilities. 

" It is proposed to determine in a series of m observations the limit of 
error J beyond which all observations involving so great an error may be 
rejected, provided there are as many as n such observations. 

" The principle upon which it is proposed to solve this problem 
is, that the proposed observations should be rejected lohen the pnAability 
of the system of errors obtained by retaining them is less than that of 
the system of errors obtained by their rejection multiplied by the proba- 
bility of making so many, and no more, abnormal observations. 



* Astronomical Joumcd (Cambridge, Mass.), Vol. II. p. 161. 
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" In determining the probability of these two systems of errors, 
it must be carefully observed that, because observations are 
rejected in the second system, the corresponding observations of 
the first system must be regarded, not as being limited to their 
actual values, but only as surpassing the limit of rejection.*' 

Let 

[i = the number of unknown quantities, 

m = the whole number of observations, 

n = the number of observations proposed 
to be rejected, 

n! =m — n, the number to be retained, 
J, J', J", . . . J^*^ = the system of errors when no observa- 
tion is rejected, 
Jj, J/, A" J . . . Jj^**'^ = the system of errors when n observa- 
tions are rejected, 
€, e^ = the mean errors of the first and second 
system, respectively, 

y =z the probability, supposed unknown, 
of such an abnormal observation that 
it is rejected on account of its magni* 
tude, 

y' = 1 — y = the probability that an ob- 
servation is not of the abnormal cha- 
racter which involves its rejection, 

X = the ratio of the required limit of error 
for the rejection of n observations to 
the mean error £, so that xe is the 
limiting error. 

The probability of an error J in the first system will be, by (14) 
and (21), 

A 1 — ^ 

^J = =: e 2€« 

ey^27r 

and the same form will be used for the second system. 

The probability of an error which exceeds the limit xe will be 
expressed by the integral (Arts. 8 and 12) 

/•a = 00 

2 I 9AdJ 
or, denoting this by if/X, 

-=^ I ^ 



2 y%A^on A* 

4x = — =: I e 2^ dJ 

ei/'2;r^A — 



'>^V74^?. 
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which, by putting t = — ^ becomes 



''=v-^fS 



and this may be found directly from Table .IX. by subtracting 

the tabular number corresponding to . = ^ from unity. 

The probability of the first system of errors, embodying the 
condition that n observations exceed the limit xe, is 



T = ipA.^A\ ipA'' 



(f(»«0/ 



1 aA«--n<c«<« 

6"" 2«« (4x)" 



c*'(27r)i*' 



in which JJ* = J» + J'* + ( J<"y ; and by (61) we have 

Sd^ = {m — [j) i?y whence 

The probability of the second system of errors is 

1 y y r 1 r 1 r 1 t^ (ItzS^^ * 






x*'(2t) 



To authorize the proposed rejection of n observations^ we 
must have 

which ^ves at once 



(^y'e*"<''-n4x)-<ry'' 



The value of y must be determined by the condition that Pj 
is a maximum, and therefore y*^''*^ y (1 — yY is a maximum. 
Taking the logarithm of this quantity, and putting its differential 
equal to zero, we obtain fpr the maximum 

n 7i' n' 
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whence 



Putting then 






^ ^ m- V (96) 



B = e*<'^-^> (4x) 



the limiting value of x, according to the above inequality, must 
be that which satisfies the equation 



(^) 






which givBs the required criterion. 

The relation of Cj to e must depend on the nature of the equa- 
tions which correspond to the rejected observations ; but it will 
give a sufficient approximation to assume that the excess of 2*4* 
over IJ^ is only equal to the sum of the squares of the errors of 
the rejected observations, which gives the equation 

(m — fi)s* — nxV = (m — fi — n) e^' 
whence 



(e, y wi — /£ — nx' 
e I m — fi — n 



m — fi 
m - 



which combined with the above equation gives 

a — nx' I T \ g^ 

=1 — - tm — n 

f^ — n \£l 
Putting, for brevity, 

'' = (e)-' (97) 



we find 



.. _ 1 = !?L:rjz:J? (1 _ ,.) (98) 



Table X.A gives the logarithms of iTand By computed by (96) 
with the aid of Table IX. We can, therefore, by successive 
approximations, find the value of x which satisfies the equations 
(97) and (98). Since B involves x, we must first assume an ap- 
proximate value of X (which the observed residuals will suggest), 
with which P will be computed by (97), and hence x by (98). 

Vol. IL— 36 
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With this first approximate value of x, a new value of log JR will 
be taken from the table, with which a second approximation to 
X will be found. Two or three approximations will usually be 
found sufficient. 

In the application of this criterion, it is to be remembered 
that it must not be used to reject n observations unless it has 
previously rejected n — 1 observations. Hence we must first de- 
termine the limiting value of x for the hypothesis of one doubtfiil 
observation, orn = 1, and if this rejects one or more observa- 
tions, we can pass to the next hypothesis, n = 2, orn = 8, &c. ; 
and so on until we arrive at the limit which excludes no more 
observations. 

The above arrangement of the tables is nearly the same as 
that given by Dr. B. A. Gould,* who was the first to prepare 
such tables and thus render the criterion available to practical 
computers. The only difference is in my table of Log. T, which 
I have found in practice to be more convenient than the corre- 
sponding one of Dr. Gould. 

Example. — "To determine the limit of rejection of one or 
two observations in the case of fifteen observations of the vertical 
semidiameters of VenuSy made by Lieut. Herndon, with the 
meridian circle at Washington, in the year 1846." In the reduc- 
tion of these observations. Prof. Pbirce assumed two unknown 
quantities, and found the following residuals {v) : 



— r.30 


— 0".24 


— 1".40 


+ 0".18 


— .44 


+ .06 


— .22 


+ .39 


+ 1 .01 


+ .68 


— .05 


+ .10 


+ .48 


— .13 


+ .20 





We have here m = 15, /4 = 2, [vv'\ = 4.2545, whence 

e» = ^— = 0.3273, e = 0".672 

We first try the hypothesis of one doubtful observation, or 
n = 1. Assuming x = 2, the successive approximations may 
be made as follows : 



* Report of the Superintendent of the U. S. Coast Suryey for 1854, Appendix, p. 
131*; also Attron. Journal, Vol. IV. p. 81. 
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2 


n 


1 


m- 


-n 


7 


m - 


p- 


- n 



n 



= 12 



Ist Approx. 

Table X.A. log T 8.404 
" « logJK 9.809 


2d Approx 
8.4044 
9.3062 


log ^ 9.095 


9.0982 


log A* 9.871 


9.8712 


log (1 — A«) 9.410 


9.4093 


I log 12 1.079 


1.0792 


log (x« — 1) 0.489 
logx« 0.610 
X 2.02 


0.4885 
0.6106 
2.020 



Hence xe = l'M6, which excludes the residual 1".40. 

We may now try the hypothesis n = 2. Commencing again 
with the assumption x = 2, we have — 



2n 




4 


m — 


n 


13 


VI — 


• At — 


n 



n 





l8t 

Approx. 


2d 
Approx. 


8d 
Approx. 


4th 
Approx. 


logT 


8.7210 


8.7210 


8.7210 


8.7210 


logi? 


9.309 


9.3622 


9.3544 


9.3558 


1 2^ 


9.412 


9.3588 


9.3666 


9.3657 


log -I' 


9.819 


9.8027 


9.8051 


9.8048 


log (1 - X') 


9.531 


9.5624 


9.5582 


9.5587 


11 , 11 

2 '**« 2 


0.740 


0.7404 


0.7404 


0.7404 


logCx"-!) 


0.271 


0.3028 


0.2986 


0.2991 


logx» 


0.457 


0.4788 


0.4755 


0.4758 


X 


1.69 


1.784 


1.729 


1.7295 



Hence xe = 0".989, which excludes the residuals 1".40 and 1".01. 
If we now try the hypothesis n = 8, we shall find, in the same 
manner, xe = 0".887, which does not exclude the residual 0".63 : 
so that the residuals 1".40 and 1".01 are in this case the only 
abnormal ones. Rejecting these residuals, we shall now find|| 
£i=0".889.* 

59. In order to facilitate the application of Pbircb's Criterion 



* For another example, in which there were four unknown quantities, and in 
which the criterion was Tery useful, see p. 207 of this Tolume. 
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in the cases most commonly occurring in practice, Table X. (first 
given by Dr. Gould) has been computed by the aid of the log T 
and log Rj according to the preceding method. 

The first page of this table is to be used when there is but 
one unknown quantity {jx = 1), or for direct observations. It 
gives, by simple inspection, the value of x' for any number of 
observations from 8 to 60, and for any number of doubtful obser« 
vations from 1 to 9. 

The second page is used in the same manner when there are 
two unknown quantities {jx = 2). 

Example. — Same as in the preceding article. — ^Having found, 
as above, e^= 0.8278, we first take from Table X. for /£ = 2 the 
value of x^ corresponding to m = 15 and n = 1, and find 

x« = 4.080, whence xV = 1.3854, xe = l'M6 

which rejects the residual 1".40. 

Then, with m = 15, n = 2, we find, from the same page, 

x« = 2.991, xV = 0.9790, xe = 0".989 

which rejects the two residuals 1".40 and 1".01. 
Passing, then, to the hypothesis n = 8, we find 

X* = 2.403, xV = 0.7865, xe = 0".887 

which does not exclude any more residuals. 

60. The above investigation of the criterion involves some 
principles, derived from the theory of probabilities, which may 
seem obscure to those not familiar with that branch of science. 
Indeed, the possibility of establishing any criterion whatever for 
the rejection of doubtful observations, by the aid of the calculus 
of probabilities, has been questioned even by so distinguished an 
astronomer as Airy.* It is easy, however, to derive an approxi- 
mate criterion for the rejection of oiie doubtful observation^ directly 
from the fundamental formula upon which the whole theory of 
the method of least squares is based. 

We have seen that the function 

'"' Remarks upon Peibcb's Criterion, Astronomieal Journal (Cambridge), Vol. IV. 
p. 137. Professor Winlock's reply to the objections of the Astronomer Royal will 
be found in the same journal, Vol. IV. p. 145. 
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"dt 



(the value of which 19 given in Table IX. A) represents, in general, 
the number of errors less than a — rt' which may be expected to 
occur in any extended series of observations when the whole 
number of observations is taken as unity, r being the probable 
error of an observation. If this be multiplied by the number of 
observations = m, we shall have the actual number of errors less 
than rt'; and hence the quantity 



..000 = 



r*[l-0O('}] 



expresses the number of errors to be expected greater than the 
limit rt'. But if this quantity is less than ^, it will follow that 
an error of the magnitude rt' will have a greater probability 
against it than for it, and may therefore be rejected. The limit 
of rejection 0/ a single doubtful ofiseruaft'on, according to this simple 
rule, is, therefore, obtained from the equation 
^ = TO [1 — eOC)] 



Qipf) = 



(99) 



If we express the limiting error under the form ze, e being the 
mean error of an observation, we shall have 

>c = — = 0.6745(' (100) 

"With the value of Q{pt') given by (99), we can find (' from Table 
IX.A, and hence Jt by (100). 

Example. — To find the limit of rejection of one of the obser- 
vations given on p. 562. "We there have m ^ 15, e — 0".572 ; 
and hence, by (99), ®{pt') = 0.96667, which ia Table IX.A cor- 
responds to ('= 3.155, whence, by (100), x = 2.128, xs ^ 1".22, 
which agrees very nearly with the limit found by Peirce's 
Criterion. 

By the successive application of this rule (with the necessary 
modifications), it may be used for the rejection of two or more 
doubtful observations, and I have, by means of it, prepared a 
table which agrees so nearly with Table X, that, for practical 
purposes, it may be regarded as identical with that table. For 
the general case, however, when there are several unfcnown 
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quantities and several doubtful observations, the modifications 
which the rule requires render it more troublesome than Peirce's 
formula, and I shall, therefore, not develop it further in this 
place. What I have given may serve the purpose of giving the 
reader greater confidence in the correctness and value of Peirce's 
Criterion. 
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■ 










n 












HI 


1 


9 


8 


4 


S 


6 


% 


8 


• 




3 


1.480 




















4 


1.912 


1. 163 


















6 


2.278 


1.439 



















6 


2!Sl6 


1.687 


1.208 












...... 




7 


I.9I0 


1.409 


1.045 














8 


3.109 


2.II2 


1.589 


1.229 

1.388 














9 


3-3»7 


2.205 
2.464 


1-753 


1.091 












lO 


3.526 


1.904 


1.531 


1.242 












11 


3.707 


2.621 


2.04c 
2.170 


1.66a 


1-373 


1.122 


••••• • 








12 


3875 


2.766 


1.78s 


1-491 
1.604 


1.2A9 
1.362 


1.018 








13 


4.029 


2.902 


2.299 


1.901 


1.145 








14 


4-173 


3.030 


2.416 


2.009 


1.709 


1.465 


1.255 


1.053 
1.163 






16 


4.309 


3.151 


2.526 


2.111 


1.807 


1. 561 


«-354 






16 


4.436 


3.264 


2.630 


2.207 


1.898 
1.985 


1.651 


1.445 


1.259 


1.080 




17 


^•5§5 


3-371 


2.719 


2.300 


1.736 
1.817 


I.|29 

1.609 


1.347 
1.428 


1.176 




18 


4.668 


3-475 


2.824 


2.389 


2.069 


1.261 




19 


4.878 


3.664 


2.914 


1-474 
2.556 


2.150 


1.895 


1.685 


1.504 
1.576 


1.341 




20 


3.001 


2.227 


1.970 


«.757 


M15 




21 


^975 


3-755 


3.084 


2.634 


2.301 


2.041 


1.827 


1.644 


1-488 




22 


5.068 


3.840 


3.164 


2.709 


1.373 


2.109 
2.176 


1.893 


1.710 


'•549 
1.612 




23 


5157 


3.923 


3.240 


2.782 


2442 


1.957 


1.773 
1.833 




24 


5.242 


4.002 


3-315 


2.852 


2.509 


2.240 


2.019 


1.671 




25 


5-3*4 


4.078 


3.387 


2.920 


*-573 


2.302 


2.079 


1.892 


1.729 




20 


5.403 


4.151 


3-456 


2.986 


2.636 


2.362 


2.137 


1.948 


'Z** 




27 


5-479 


4.222 


3-588 
3.651 


3.049 


2.697 


2.420 


1-194 


2.001 
2.056 


1.838 




28 


5-$S» 
5.622 


4.291 


3. Ill 


2.756 


1-477 


1-149 


1.891 




29 


4.358 


3.171 


2.813 


2.586 


2.302 


2.108 


1.941 




30 


5.690 


4.422 


3.712 


3.229 


2.869 


2-354 


2.158 


1.990 




31 


5756 


4.484 


3.772 


3.285 


2.921 
2.976 


2.618 
2.68q 
2.718 
2.787 

1.834 


2.404 


2.207 


2.038 
2.085 




32 


5.820 


4.545 
4.604 


3-5;9 


3-340 


2.454 


2.255 




33 


5.882 


3.884 


3-394 
3-446 


3.028 


2.502 


2.303 


2.130 




34 


S-94» 
o.ooi 


4.661 


3.939 


3.078 


2.549 


1.347 


1.175 

2.21 1 




35 


4.717 


3.992 


3-497 


3.127 


2.594 


2.392 




30 


6.058 


4-771 


4.044 


3-547 


3-174 


2.880 


4.619 
2.681 
2.726 


2.436 


2.261 




37 


6.II1 


4.823 


4.095 


3-595 
3.689 


3.221 


2.926 


2.478 


2.302 




38 


6.167 


4.874 


4-144 


3.267 


2.970 


2.520 


*-343 
2.383 




30 


6.219 


4-9*5 


4.192 


3.312 


3.013 


2.768 


2.c6i 
2.601 




40 


6.270 


4-974 


4-»39 


3-734 


3.356 


3055 


2.809 


2422 




41 


6.320 


5.022 


4.285 


3-779 


3-398 


3.097 
3.138 


2.888 


2.640 


2.460 




42 


6.369 


5.069 


4.331 


3.822 


3-440 
3-481 


2.678 


a-497 




43 


6416 


5.114 


4.375 


3.86J 
3.906 


3.178 


2.927 


2.716 


a-534 




44 


6.461 
6.508 


5-«59 


4.418 


3.521 


3.217 


2.965 


1-753 
2.789 


xm 




45 


5.202 


4.460 


3.947 


3.561 


3.255 


3.002 




40 


6.552 


5-H5 
5.287 


4.501 


3-987 


3.600 


3-*93 


3-039 


*!;5 


2.641 




47 


6.596 


4-54* 
4.581 

4.620 


4.026 


3.638 


3.366 


3-075 


2.860 


a.67j 
2.708 




48 


5-529 


^•^ll 


4.065 


3675 


3.110 


2.89J 
2.928 




1 40 


6.681 


5.368 


4.103 


3.712 


3.401 


3-145 


2.741 




50 


6.720 


5.408 


4.657 


4.140 


3.748 


3436 


3179 


2.962 


*.774 




51 


6.761 


5-447 
5-4*4 


4-695 


4-176 


3.784 


3-471 


3.211 
3.246 


a.994 


2.806 




52 


6.800 


4.768 


4.212 


3.819 


3-50J 
3-538 


3.027 


2.838 
2.869 




53 


6.818 


5.522 


4-*47 
4-»82 


3?53 


3179 


3059 




54 


6.876 


5-559 


4.804 


3-887 


3-571 
3.603 


3.311 


3.090 


2.899 




55 


6.913 


5-595 


4-839 


4-316 


3.920 


3-342 


3.121 


2.929 




56 


i'Ht 


5.610 
5.665 


4.873 


4-349 
4.382 


3-951 


lUl 


3.373 


3.151 
3.181 


^^2! 




57 


6.986 


4.907 


3.984 
4.016 


3.404 




58 


7.021 


5.699 


4.941 


4-415 


3.697 


3-434 
3.463 


3.210 


3.017 




50 


7.056 


5.766 


4.974 
5.006 


4-447 


4.047 


3-718 


PI 


3.046 


}/ '^l 


7.090 


4-478 


4.078 


3.758 


3.492 


3-074 



5^« 



' 
















1 


^^^H 








TABLE X. Peirce'a 


Criterion. 




H 








Values or 


«'F0B^=2. 






^ 


f 




















I 




1 


a 


8 


4 


5 


e 


1 


B 


« 




! 4 


..484 


















ft 




« 


»35 






















' 6 

T 


1.150 


1 


479 


;:;;« 



















8 










...63 






































10 


J.I4» 


a 


»77 


i-77> 


1.439 


1.191 














11 


3*17 


1 


440 


1.911 


..566 


1.310 


l.ogf 





































3.78; 






1.171 






1.310 


















1.19c 


1.91c 




1.40! 




1.045 










1 


991 


1,403 




1.7x7 


..501 


i.jii 








16 


4.1. s 


1 


109 


1.510 




1.819 


1.589 


\X 


1.119 


1.070 






tir. 










1.90; 












1 






1.195 


1.991 


'■75! 


1.557 


































4.«U 


3 


5ii 


5.89= 


1.465 


1.150 




1.703 


1.531 


..380 




31 






i.q 








.,976 


1.771 


1.598 








4-Hs 










i.igfi 




..83! 


































S..6, 




(74 


i.iii 


zAl 






1.964 


1,785 


•h' 






S"SS 


3 


953 


J.ii6 


i-hs 


1.501 


i.ljr 


1.014 


1.84; 






26 


S."!« 


4.D19 


3-357 


1.901 


1-565 


1.199 


1.0S1 


1.900 


'■743 




























«?J 


*-'7-f 


1'49» 






1.415 


t-llt 




























5.559 


4-J09 


J.S'9 


3.150 




1.515 


1.300 








31 






J.6S0 


3.10B 


1.856 


I.57S 




1,160 


1.996 




33 






J'739 
3.796 


3.164 




dy 


1.401 


i.iot 








5.740 




3-319 






1.445 


1,155 










4-i5S 






3,01 = 


1,715 




1.301 


1,131 




35 


5.864 


4.613 


J.JOb 


3.414 




»-777 


»-S4: 


»-345 






3S 




4.669 








1.81* 


i.j8( 


1.389 


1.119 








*-7»3 




3'5i3 








1431 






38 












1,9.4 

1.95? 




1.474 






30 










3.150 


1.717 


1.51; 


1,341 




411 


6.145 




4-159 


1-665 


3-»94 


3.001 


1-759 


*-5SS 




] 


41 






4.106 




V"8 




1.800 








6.147 


4.975 


4.151 


3.755 




j.o«4 






1-457 


43 










J.41J 












44 


*'34S 






lit: 






1,917 


1.705 


1.530 




45 


6.391 


5.113 




3-S°3 


3-"] 


»-95S 








46 


6.4ja 


S''S7 


4.416 


3.913 


3543 


3.141 


1,991 


1.781 


1.601 




47 






















48 


6.517 


\'.r, 


4.50S 


4.1101 


j.^.S 




3.064 


Hn 








6.S7C 








3-35' 


3-099 




1.703 




50 




S'3»3 










3.34 


1.910 






51 


6,653 


5,361 


4.616 




3.718 


3.411 


,..68 


1.953 


1.76B 




53 




5.40. 




4.1 <■ 


SS) 




3.10. 








53 




S.+4<: 








3-134 




1.831 




54 






















55 




S-SI5 




4-«57 


3-«67 


3SSi 




j,o5o 


1,891 




56 


6.S48 
6,lti5 


S55' 


4.807 


4.191 


3.900 


3.588 


3.319 


3.111 


1-911 




57 




















58 








4.357 






3.39° 


3.171 






5& 


6-951 




4 9°! 


♦■39= 






3.4'9 


3.100 


3.009 




1 flU 


6-993 




4.941 


4,411 


4.017 


3711 


3.448 


3.119 


3.037 














591 










k 


L 





















TABLE X. A. Feiroe'a Criterion. 

Log T. 













« 






























1 


i 


8 


t 


4 


6 


I 


8 


» 


J 


9-397i 








1 




!l 


3 


958S3 
9.3979 
9.169J 


















9-6744 

9.5119 

9.1338 












5 


B.9'14 
8.8159 


9.7183 

9.1693 










s 


9.7651 
9.6381 

9+631 

9-3979 












ill; 












8 




9-5853 
9,5119 


g-Bijo 









8.9134 


9.70+1 
9.61.0 


98196 

9.7153 




10 


9-1157 


9-8+3' 


11 


8.5+47 


g.867S 


9.0669 


9.1171 


9-3+' 7 


9-4i'4 


9-SS17 


9.6|oi 


9.7+83 


12 


B.joji 


i^^V> 


9,0131 


9.1707 


9-1911 


9-3979 


9.49+3 


9.5853 


967++ 
9,6ii3 


13 


8.46^9 


M- 


S-983+ 


9.1188 


9.1+77 


9.3506 


9.4433 


9.519a 


14 


8.43S! 


8-7S3> 


8.9470 


9.0906 


9.1074 


9.J080 


9.3979 


9.4810 


9-5597 


IS 


a.4044 


8.7110 


8-913+ 


9-OSS5 


9.. 707 


9.169, 


9.3570 


9+37+ 


9.5.19 


l« 


!■"!* 


8.69.0 


g.88ii 


9,013. 


9..368 


9.1338 


9-3197 


9-3979 


9.47.0 


n 


B.)4Sj 


8.6619 


8.8;ji 


8-9?3° 


9.1055 


910.. 


9-185+ 


9.3619 


9.4318 


18 


!-^"JZ 


8.6-165 


8.8159 


8.9^+g 


9.0761 


9-1707 


9-1537 


9.3187 


93979 


19 


3.19S6 


8.6117 


8.8003 


8-9383 


9.0489 


^.1413 


9.11+1 


9.1980 


9.365* 


20 


8.1757 


8.58.81 


8.776. 


8.913+ 


9.0131 


9-' "57 


9.1966 


9.1693 


9-3359 


21 


8.1540 


g.S6s9 


8.7531 


8,8898 


8.9988 




9.1707 


9.141+ 


9.3080 


2S 


8.2,3, 


8-5+47 


8.7315 


Ifl^ 


8.9758 


9:0669 


9-1463 


9.1.71 


9.18.8 


33 


8.1156 


g.ji+i 


a.7107 


8,8+61 


8.95+0 


9-0+4S 


9.1131 


9.1933 


9.1571 


24 


8.1947 
8.1766 


8.5051 


8.6910 


8,8159 


8-9 3 31 


9.0131 


9.1011 


9.1707 


9-1338 1 


SS 


8.48^7 


8.6711 


8,8066 


8.9134 


9-ooiB 


9.0803 


9.1+91 


9,1117 ] 


2« 


8.1591 
8,I4;S 


8.4689 
8.4519 


8.65,9 
8.6365 


8.788. 
8.7703 


a?K 


\tli 


9,060+ 
9.0+1+ 


9.1091 


9.1907 
91707 


28 


8. 1 264 


8-435+ 


B60'' 


l-^ni 


8.8588 


8.947° 


9-0131 


9.0906 


9..S.6 


20 


g.M.39 


8.4197 




a. 7368 


8.8+1. 


8.9199 


g-°2Ji 


9.0717 


9,1331 


30 


8.0959 




8,5881 


8.7110 


8.8159 


8.9.34 


S.ggSS 


9-0555 


9."57 


31 


8-08.4 


8,3897 


8.5731 


8-7057 




8-8975 


8,9716 


9.0390 


g.DBii 


32 


8.0674 


8.375+ 


a-5SS7 


8.691a 


8-7954 


8.8811 


8.957. 


9,0131 


9,oii6 1 


3J 


8.o;,g 


8.3617 


8-S++7 


8.6767 


8,7805 


8,8675 


8,9410 


9,0078 


9.0669 


34 




8.3483 


8,53.1 


8.6619 


8.76GI 


8.8531 


8,9175 


8-9930 


9,0518 ■ 


35 


8"oi79 


8-3353 


8.5179 


8-6495 


8.75J1 


8.8393 


8.9134 


8.97S6 


9.0371 


38 


a.0155 


8,3117 


8,505. 


8,6,65 


S7400 


8,8159 


8.R998 


8.96+8 


9,0131 : 


37 


8.0D34 


8.3105 


8.4917 


8-6139 


B.7171 


8,8.19 


8.8665 


8,95.3 


9,0095 1 


38 


7-9917 


81986 




8.6,17 


8.7 48 


8,8003 


8.8737 


8.9 83 


K.9961 


30 


7.98QJ 


8,^70 


g:j6g9 


8.5998 


B.7017 


8.7881 


8.8613 


B.9157 


89834 


40 


,.9691 


8.1757 


8,4575 


8.5881 


8.6910 


8,7761 


B.8491 


8.9,3+ 


8,970:; 


41 


7.9 5B3 


8.1647 


8.4463 


8.5769 


8.6795 


8.764; 


8.837+ 


l%t\t 


B,95i^S 


ta 


7-9+77 




8-+35J 


8.5659 


8.668+ 


8.753* 


B.ai;, 


8-9-J7C 


43 


7-957J 


8.1435 


8.41+9 


8.5551 


8-6575 


8.7411 


8.8 148 


8,8785 




44 


7-9^7* 


8.1133 


|.+'+S 


8-54+7 


a.6+69 


B.7315 


B.8039 


8.867; 


tli\\ 


43 


7-9 '74 


8.1133 




3-53+5 


S.6365 


B.7110 


«-7933 


8.8567 


8-9- 34 


46 


7-9°77 


a.1.36 


8-394! 


8.51+5 


8.6164 


8,7.07 


8.7819 


8.B461 


89C1- 


47 


Hf-i 


8.1040 


8,38+9 


8,51+7 


a.6165 


8,7007 


a.7716 


8,8360 


IJ914 


48 


7,8890 


8.1947 


8-3754 


a.505. 


8.6069 




8.7619 


8,8159 




40 


7.8800 


8.1I55 


8.3661 


8.4958 


a. 5974 




8.7531 


8,8.61 


ais;.'^ 


H) 


7.S71. 


8.1766 


8.3571 


8.1^^7 


t\llX 


8.671. 


8.7+38 


8,8066 


8,S6i(. 


51 


7.8614 


8.1678 


8.3+83 


^■+777 


8.5791 


8,6619 


8-73+5 


\Wu 


8,8531 


52 


7-8339 


8.1591 


8-3396 


8.4689 


8.5703 


8-6539 


8,7154 


8,8439 


53 


7.8+55 


8.,;o8 




8.4603 


8,?6i6 


8.6451 


8.716? 


8-779" 


8,83+^ 


54 


7-8374 


8.1+15 


8'3"7 


8,4519 


8-553° 


8.6365 


8.7079 


8-7703 


8,8359 


5S 


7-8193 


8-13+4 


8,31+5 


B-J+3S 


8-5447 


8.618? 


8.69^1 


8.76,7 




56 


7.8X.4 


a. 1 164 

8.1 i8g 


8,3065 


8.+35; 


8,5365 


8.6198 


8.69.0 


g.7531 


8,8087 1 


57 


7.8157 


8.1986 


8-+17; 


8,5184 




8.6818 


V,% 


8,8003 1 


58 


7.8060 


8.1109 


8,1503 


8.+ . 97 


B-510J 


b:6o37 


llltl 




50 


7.7986 


8-1033 


8.Jji 






8-5959 


■ .7.J8 


87*|0 


OU 




8-0959 


8.1757 


g.+o++ 


Sijos. 


8.58^1 


8.6590 


8.7110 



^ 




TABLE X. A. Feirce'a Griterioii. _ ^^| 








Log T. 


■ 






■ 


1 


























1 


i 


3 


4 


5 


e 




J 


8 


9 




ai 


I'.jytt 


8.08 !6 


8.2684 


1.3970 


8-4977 


8.5806 




6514 


a-7.33 


8.7684 




vi 


8.0814 


8 






3S97 


8.4903 


8.5731 




6439 


8-7057 


8.7607 








7.769a 




8 


1540 


( 


3S1S 


8.4630 


8.5659 




6365 


8. 698 3 


8.7531 






64 


7.7619 


8.067. 


8 


1470 


8 


3754 


8-4759 


8.5587 




6293 


B.6gio 


8.7458 






69 


7.7561 


3.0606 


8 




8 


3685 


S.4689 


8.5516 




6111 


8.6838 


e.7386 






68 


7-7495 


80538 


8 


»333 


8 


36,7 


8.4610 


8.5447 




6.51 


8.6767 


8.73.5 






6T 


7.7^19 


B.0+72 




1266 


! 


3549 


S.4551 


8.5378 




60E1 


8,6697 


8.7144 






68 


7-7364 


S.0407 


8 






3483 


S.4485 


8-S3" 
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